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Abstract 

Cosmological observations suggest that the Universe is undergoing an accelerated phase of 
expansion driven by an unknown form of matter called dark energy. In the minimal stan- 
dard cosmological model that best fits the observational data the dark energy is provided 
by a cosmological constant term. However there is currently no convincing theoretical 
explanation for the origin and the nature of such an exotic component. In an attempt 
to justify the existence of the dark energy within the framework of particle physics the- 
ories, several scenarios have been considered. In this thesis we present and discuss the 
phenomenological aspects of some of these dark energy models. We start by reviewing 
the cosmological measurements that give direct and indirect evidence for the dark energy. 
Then we focus on a class of theoretical models where the role of the dark energy is played 
by a minimally coupled scalar field called quintessence. For the sake of simplicity we place 
special emphasis on two of these models, the Inverse Power Law times an Exponential 
potential and the Two Exponential potential. We then consider the effect of scalar field 
fluctuations on the structure formation process in these two models. By making use of 
the cosmological distance measurements such as the supernova luminosity distance and 
the position of the acoustic peaks in the CMB power spectrum, we constrain the shape 
of a general parameterized quintessence potential. We find that by the present time the 
scalar field is evolving in a very flat region or close to a minimum of its potential. In 
such a situation we are still unable to distinguish between a dynamical model of dark 
energy and the cosmological constant scenario. Going beyond constraining specific classes 
of models, we develop a model independent approach that allows us to determine the 
physical properties of the dark energy without the need to refer to a particular model. 
We introduce a parameterization of the dark energy equation of state that can account 



for most of the proposed quintessence models and for more general cases as well. Then we 
study the imprint that dark energy leaves on the CMB anisotropy power spectrum. We 
find that dynamical models of dark energy produce a distinctive signature by means of the 
integrated Sachs- Wolfe effect. However only models characterized by a rapid transition 
of their equation of state can most likely be distinguished from the cosmological constant 
case. By using a formalism to model localized non-Gaussian CMB anisotropics, we com- 
pute analytical formulae for the spectrum and the bispectrum. The use of these formulae 
in specific cases such as the SZ signature of clusters of galaxies provide an alternative 
cosmological test. 
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Introduction 



The set of astrophysical observations collected in the past decades and the theoretical and 
experimental developments in high energy physics have provided the natural framework 
that defines Cosmology as a scientific discipline. The identification of the 'Hot Big-Bang' 
scenario as a paradigm has been of crucial importance for the beginning of Cosmology 
as a modern science. In fact it has allowed us to address a number of questions about 
the nature and the evolution of the Universe that otherwise would have been the subject 
of investigation of philosophers and theologians. Since this paradigm has been accepted 
by the majority of the scientific community, more specific and detailed studies have been 
undertaken in order to extend the validity of the paradigm to a wider class of phenomena 
such as the formation of the structures we observe in the Universe. As result of this intense 
activity, that the philosopher of science T.S. Kuhn would define as normal science inves- 
tigation [1], the initial paradigm of Cosmology has been extended in order to include the 
inflationary mechanism and dark matter, two necessary ingredients to explain a number of 
issues that arise within the 'Hot Big-Bang' scenario. This extended paradigm is extremely 
successful and recent measurements in observational cosmology have widely confirmed its 
prediction. It is very remarkable that long before the recent developments of the cosmo- 
logical science, T.S. Kuhn identified the basic steps that a scientific discipline follows in its 
evolution, steps that applies to modern Cosmology too. In particular he has pointed out 
that there are phenomena which evading an explanation of the paradigm are the subject 
of 'extraordinary' investigations, in opposition to the 'ordinary' normal science activity. 
Most of the time these studies lead to a crisis of the underlying paradigm and trigger 
what he has called a 'scientific revolution'. In this light the discovery that the Universe is 
dominated by a dark energy component that accounts for 70% of the total matter budget 
belongs to this class of phenomena. This is the subject of this thesis. In Chapter ^ we will 
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describe the observational evidence of dark energy and in Chapter|21we will review some of 
proposed explanations. Prom Kuhn's point of view these explanations would represent the 
attempt to force Nature to fit within the 'Hot Big-Bang' paradigm. In fact these models 
fail to succeed since they manifest inconsistencies with the expectations of particle physics 
theories that are included in the standard cosmological paradigm. The solution to this dif- 
ficulty most probably will need a new paradigm. This necessity will be more urgent if the 
observational data will indicate a time dependence of the dark energy properties. In this 
perspective the aim of this thesis is to investigate some of the phenomenological aspects of 
minimally coupled quintessence scalar field scenarios. In Chapter |2l we describe the evo- 
lution of fluctuations in the quintessence field. In Chapter |1] we discuss the constraints on 
a general class of quintessence potentials obtained from the analysis of the position of the 
acoustic peaks in the Cosmic Microwave Background (CMB) anisotropy power spectrum 
and the Sn la data. In Chapter [S] we describe some of the methods used to constrain the 
dark energy. We then introduce a model independent approach that allows us to study 
the full impact a general dark energy fluid has had in Cosmology. In fact this fluid de- 
scription, based on a very general parameterization of the dark energy equation of state, 
allows us to infer the dark energy properties from cosmological observations, instead of 
constraining specific classes of dark energy models. In Chapter El we study the effects dark 
energy produces in the CMB anisotropy power spectrum and show that clustering of dark 
energy leaves a distinguishable signature only for a specific class of models. In Chapter [7] 
we present the results of preliminary work that aims to develop alternative cosmological 
tests using the non-gaussianity produced by localized sources of CMB anisotropics. We 
hope that the work reviewed in this thesis will provide the basis for those 'extraordinary' 
investigations that will help us in developing the new paradigm that modern Cosmology 
needs. 



Chapter 1 

The dark side of the Universe 



The recent results obtained in different areas of observational cosmology provide an as- 
tonishing picture about the present matter content of the Universe. The accurate mea- 
surements of the Cosmic Microwave Background radiation (CMB) give strong evidence 
that the curvature of the space-time is nearly flat. On the other hand the analysis of 
large scale structure surveys shows that the amount of clustered matter in baryonic and 
non-baryonic form can account only for thirty per cent of the critical energy density of the 
Universe. In order to be consistent these two independent analyses require the existence 
of an exotic form of matter, that we call dark energy. More direct evidence is provided 
by the Hubble diagram of type la supernova (Sn la) at high redshifts. It suggests that 
the Universe is undergoing an accelerated expansion sourced by this dark energy compo- 
nent that is characterized by a negative value of its equation of state. A lot of criticism 
has been levelled to these measurements. In fact the physics of Sn la is still a matter of 
debate and consequently their use as standard candles has not yet convinced the whole 
astronomical community. In spite of such an important issue it is worth underlining that 
the combination of CMB and supernova data constrains the amount of clustered matter 
in the Universe to a value that is consistent with the large scale structure observations. As 
we shall review in this chapter, indirect evidence for a dominant dark energy contribution 
comes mainly from the combination of CMB results and the constraints on the density 
of baryons and cold dark matter. The simplest explanation for the dark energy would 
be the presence for all time of a cosmological constant term A in Einstein's equations of 
General Relativity. However evidence for a non vanishing value of A raises a fundamental 
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problem for theoretical physics, (for a review of the subject [2-7]). In fact, as we shall 
see later, it is rather difficult to explain the small observed value of A from the particle 
physics point of view. In this chapter we will briefly review some historical developments 
of the dark energy problem. We will introduce the standard cosmological model and the 
equations describing the expansion of the Universe. Then we will discuss the build up of 
observational evidence for dark energy. 

1.1 A historical introduction 

The cosmological constant was initially introduced by Einstein in the equations of General 
Relativity (GR) as a term that could provide static cosmological solutions [8]. Motivated 
by the observed low velocities of the stars, he assumed that the large scale structure of the 
Universe is static. We should remind ourselves that the notion of the existence of other 
galaxies had been established only a few years later. Besides, Einstein believed that the GR 
equations had to be compatible with the Mach's principle, which in a few words states that 
the metric of the space-time is uniquely fixed by the energy momentum tensor describing 
the matter in the Universe. For this reason he thought the Universe had to be closed. 
These two assumptions were, however, not compatible with the original form of the GR 
equations. In fact a matter dominated closed Universe is not static, therefore he needed 
to introduce a term leading to a repulsive force that could counterbalance the gravity. In 
such a case he found a static and closed solution that preserved Mach's principle. But 
in the same year, 1917, de Sitter discovered an apparent static solution that incorporated 
the cosmological constant but contained no matter [9]. As pointed out by Weyl, it was 
an anti-Machian model with an interesting feature: test bodies are not at rest and an 
emitting source would manifest a linear redshift distance relation. Such an argument 
was used by Eddington to interpret Slipher's observations of the redshift of spiral nebula 
(galaxies). Subsequently expanding matter dominated cosmological solutions without a 
cosmological constant were found by Friedman [10,11] and Bubble's discovery of the linear 
redshift distance relation [12] made these models the standard cosmological framework. 
The cosmological constant was then abandoned. For some time a non vanishing A was 
proposed to solve an ^age problem'. Eddington pointed out the Bubble time scale obtained 
by using the measured Bubble constant was only 2 billion years in contrast with the 
estimated age of the Earth, stars and stellar systems. In the 1950s the revised values of 
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the Hubble parameter and the improved constraints on the age of stehar objects resolved 
the controversy and the cosmological constant became unnecessary. However in 1967 it 
was again invoked to explain a peak in the number count of quasars at redshift z = 2. It 
was argued that the quasars were born during a hesitation era , at the transition between 
the matter and the A dominated era. More observational data confirmed the existence of 
this peak and allowed for a correct interpretation as simply an evolutionary effect of active 
galactic nuclei, with no necessity for A. We shall discuss the theoretical implications of 
the cosmological constant in Chapter 2, here we would like to stress that in the past few 
decades the A term has played the role of a fitting parameter necessary to reconcile theory 
and observations which has been discarded every time systematic effects were considered. 
It is therefore natural to ask the question if today we are facing a similar situation. In 
what follows we will try to show that this turns not to be the case and the dark energy is 
indeed most likely to be present in our Universe. 

1.2 Standard cosmology 

The Einstein field equations are: 



where K^y is the Ricci tensor, R is the Ricci scalar, A is the cosmological constant term 
and T^y is the matter energy momentum tensor which determines the dynamics of the 
Universe. When different non interacting sources are present the energy momentum tensor 
is the sum of the energy momentum tensor of each of the sources. Assuming an isotropic 
and homogeneous space-time the large scale geometry can be described by the Priedman- 
Robertoson- Walker metric: 



where a(i) is a function of time (called the scale factor) and k = 0, ±1 sets a flat, open 
(-1) or close (+1) geometry. Spatial homogeneity and isotropy implies that the energy 
momentum tensor of each component is diagonal: 




(1.1) 




(1.2) 




diag{pi{t),pi{t),pi{t),pi{t)) 



(1.3) 
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where Pi{t) is the energy density and Pi{t) is the pressure of the i-th. matter component 
(radiation, baryons, cold dark matter, etc.). In the FRW metric the Einstein equations 
with a mixture of different matter components are the Friedman equations: 



tt2 I 8ttG st-^ a k . ^. 



a J 3 ^-^ 3 a 
- = -_^(p„+3,.)+3^ (1.5) 

i 

We define the density parameters Q.i = Pi/pc, = ^/'^Pc and = —k/H'^o? where 
Pc = 3H'^ /SnG is the critical energy density. Then the Friedman equation Eq. (|1.4j) can 
be rewritten as: 

i-nk = Y.^i = ^tot, (1.6) 

i 

showing that the spatial curvature is fixed by the total matter content. Since the different 
components do not interact with each other, their energy momentum tensor must satisfy 
the energy conservation equation T^.^, = 0. Hence in addition to the Friedman equations 
the evolution of the energy density of each matter component is given by: 

Pi = -3H{p^+p,). (1.7) 

It is worth remarking that, as has been stressed by T. Padmanabhan [16], 'absolutely no 
progress in cosmology can be made until a relationship between pi and pi is provided in 
the form of the functions Wi{a) '. In fact once these relations are known, we can solve 
the dynamical equations and make predictions about the evolution of the Universe, that 
can be tested by cosmological observations. If the matter components consist of normal 
laboratory matter, then the knowledge of how the matter equation of state w evolves 
at different energy scales is provided by particle physics. At present the behaviour of 
matter has been tested up to about 100 GeV, in this domain the relation between energy 
density and pressure can be taken to be that of an ideal fluid, pi = WiPi, with w = 
for non relativistic matter and w = 1/3 for relativistic matter and radiation. However 
if a cosmological model based on conventional matter components fails to account for 
cosmological observations, we could interpret this fact as a failure of the cosmological 
model or as a signal for the existence of a source not seen in laboratories. For instance 
the cosmological constant term behaves as a perfect fluid with negative pressure. This 
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can be seen rewriting the A term in Eq. H1.4p and Eq. 1)1. 5|) as an energy density and a 
pressure term. Then one finds p\ = —p\ = — A/(87rG). The effect of such a component 
on the expansion of the Universe can be seen from Eq. ()1.5)) , and the value of deceleration 
parameter today is 



and we have neglected the radiation. For VL^ > 0^/2 the expansion of the Universe is 
accelerated since qq < 0. Hence in a Universe dominated by the cosmological constant 
the expansion is eternally accelerating. In summary the dynamics of our Universe is ob- 
servationally determined by two geometrical quantities, the Hubble parameter Hq, which 
provides us with a measure of the observable size of the Universe and its age, and the de- 
celeration parameter qo which probes the equation of state of matter and the cosmological 
density parameter. 

1.3 Observational evidence 

Different cosmological tests can be used to constrain the geometry and the matter content 
of the Universe. We shall briefly review the latest limits on and ^^A obtained by recent 
experiments in cosmology. 

1.3.1 CMB anisotropies 

During the last few years an avalanche of balloon and ground experiments, together with 
the most recent WMAP satellite observatory have measured the small angular temper- 
ature fluctuations of the Cosmic Microwave Background Radiation. Such measurements 
have detected a series of acoustic peaks in the anisotropy power spectrum and confirmed 
early predictions about the evolution of pressure waves in the primordial photon-baryon 
plasma [17, 18]. The specific features of such peaks are sensitive to the value of the cosmo- 
logical parameters, in particular to ^Itot, and the scalar spectral index n. The sensitivity 
to the curvature of the Universe however does not allow us to constrain independently 
and r^A) that are consequently degenerates. The earlier analysis of the Boomerang ex- 
periment [19-21] found rjfc ~ and the latest data released constrain the total energy 
density to be Qtot = ^■^'^^o'ol t^^l- Such a result is consistent with the ones found 
by other CMB experiments. For instance, the data from MAXIMA-1, another balloon 
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experiment, when combined witii the COBE-DMR data suggest Qtot = l-OOioio I^"^]- 
Similarly the two ground experiments, DASI and CBI provide i^tot = 1-04 it 0.06 [24] 
and Qtot = 0.99 ± 0.12 [25] respectively. Recently three more groups, ARCHEOPS [26], 
VSA [27] and ACBAR [28] have released their data finding similar results. The constraints 
on the baryon density are in good agreement with the prediction of the Big-Bang Nucle- 
osynthesis (BBN) and the scalar spectral index is found to be of order unity, as predicted 
by generic inflationary paradigms. However CMB alone poorly determines Q\ and a van- 
ishing cosmological constant cannot be excluded at 2a. Nonetheless due to the strong 
constraint on the curvature of the Universe, it is reasonable to restrict the data analysis 
to the flat cosmological models {^tot = !)• In this case, assuming the so called 'HST 
prior' on the value of the Hubble constant, h = 0.71 it 0.076 [29], then all the CMB data 
constrain the cosmological constant density parameter to be f^A = 0.69ibQ;Qg [25]. The 
WMAP satellite provided CMB data with such an high level of accuracy that is worth 
mentioning a part. The experiment has measured CMB anisotropics in different frequency 
bands, allowing for an efficient removal of the foreground emissions. The measurements 
mapped the full sky in the unpolarized and polarized components providing an accurate 
determination of the temperature power spectrum (TT) and the temperature-polarization 
cross-correlation spectrum (TE) [30]. The position of the first peak in the TT spectrum 
constrain the curvature to be Qk = 0-030ibg;Q25 [31]. The combination of WMAP data 
with ACBAR and CBI, 2dF measurements and Lyman a forest data find the best fit cos- 
mological parameters: h = 0.71ibQg|, the baryon density Ofo/i^ = 0.0224 it 0.0009, the dark 
matter density ^mh'^ = 0-0135±Q;QQg, the optical depth r = 0.17 ±0.04, the scalar spectral 
index n = 0.93ib0.03 and the amplitude of the fluctuations as = 0.84ib0.04 [32]. The value 
of T comes from an excess of power on the large angular scales of the TE spectrum [33]. 
This signal cannot be explained by systematic effects or foreground emissions and has a 
natural interpretation as the signature of early reionization, most probably occurred at 
redshift z w 20. This conflicts with the measurements of the Gunn-Peterson absorption 
trough in spectra, which indicate the presence of neutral hydrogen at redshift z ~ 6 [34] . 
Therefore we have evidence for a complex ionization history of the Universe, which most 
probably underwent two reionization phases, an early and a late one. Of particular inter- 
est is the running of the scalar spectral index that provides a better fit to the data when 
WMAP is combined with small angular scale measurements such as ACBAR, CBI, 2dF 
galaxy survey and Lyman a. Another interesting finding of the WMAP TT spectrum is 
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the lack of power at low multipole. In particular the quadrupole and the octupole are 
suppressed compared to the expectation of the best fit ACDM model. It has been claimed 
that such suppression could be the signature of new physics [35,36]. 

1.3.2 Clustering of matter 

The cosmological structures we observe today have been formed by the gravitational am- 
plification of small density perturbations. The amount of such inhomogeneities at different 
cosmological scales is measured by the matter power spectrum. This is estimated from 
the statistical analysis of a large sample of galaxies and provides a measurement of the 
amount of clustered matter in the Universe. Recently two large galaxy surveys, the 2dF 
Galaxy Redshift Survey [37] and the Sloan Digital Sky Survey [38], have probed inter- 
mediate scales (10 — 100 Mpc). The fit to the power spectrum data of the 2dF yields 
= 0.20 lb 0.03 and the baryon fraction Qf,/il.rn = 0.15 it 0.07 [37]. Such a low value 
of gives indirect evidence for a large non vanishing cosmological constant contribution 
when this LSS data is combined with the CMB. A joint analysis of the CMB and 2dF 
data indicates 0.65 < Q\ < 0.85 at 2a [39]. An independent estimate of is provided 
by the peculiar velocities of galaxies. In fact mass density fluctuations cause galaxy mo- 
tion with respect to the Hubble flow. Such a motion reflect the matter distribution and 
therefore is sensitive to Qm- The analysis of the Mark III and SFI catalogs constrain 
Qrn = 0.3 ± 0.06 [40]. Low values of Q m are also indicated by the studies of cluster of 
galaxies, where it is assumed the amount of matter in rich clusters provides a fair sam- 
ple of the matter content of the Universe. Recent surveys have precisely determined the 
local X-ray luminosity function. Using the observed mass-luminosity relation, the cluster 
mass function has been compared with the prediction from numerical simulations. This 
analysis constrains Qm. < 0.36 at la [41] (see also [42]). This result is in agreement with 
the limits found by an alternative study from which, 0.1 < 0,m < 0.5 at 2a [43]. Another 
way of estimating the amount of dark matter is to measure the baryon fraction fb from 
X-ray cluster observations. In fact the ratio of the baryonic to total mass in cluster should 
closely match the ratio Qfj/^m- Therefore a measurement of fb combined with accurate 
determination of Qb from BBN calculation can be used to determine Qm- Using such a 
method it was found ~ 0.32 for h 0.7 [44]. This is in agreement also with the 
value obtained by a study of the redshift dependence of the baryon fraction, that indicates 
= 0-3±o'o3 [45]. Similarly a different analysis based on gravitational lens statistics 
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provides Qm = O.SlzbgJ^ [46]. 
1.3.3 Age of the Universe 

The Friedman equation Eq. 1)1. 4|1 can be integrated to obtain the age of a given cosmological 
model: 

da , , 



/o ay^ (1 - Qrn - ^A)/aP' + ^m/a^ + ^^A 
The numerical solutions are shown in figure ITTTl The solid lines correspond to increasing 
value of HqIq in the $7^ — Vt\ plane. It is easy to see that for fixed values of $7^ the age 
of the Universe increases for larger values of f^A. Since matter dominated cosmological 
models are younger than globular clusters, in the past few decades the possibility of an 
'age problem' has been a matter of debate. The presence of a cosmological constant can 
alleviate such a problem. However a key role is played by the Hubble parameter, in fact 
low values of Hq increase to- The age of globular clusters is estimated to be about 11. 5 it 1.5 
Gyr [47], therefore a purely matter dominated Universe (fim = 1 and f^A = 0) cannot be 
excluded if /i < 0.54. However such a low value of h seems to be inconsistent with the 
accepted value of /i ~ 0.71 it 0.07 [29]. It is worth mentioning that the determination 
of the age of high redshift objects can be used to constrain Oa by studying the redshift 
evolution of the age of the Universe [48] . 

1.3.4 Supernovae la and luminosity distance measurements 

Supernovae type la are violent stellar explosions, their luminosity at the peak becomes 
comparable with the luminosity of the whole hosting galaxy. For such a reason they are 
visible to cosmic distances. These supernovae appear to be standard candles and therefore 
are used to measure cosmological distances by means of the magnitude-redshift relation: 

m-M = 51ogio^ + 25, (1.10) 
Mpc 

where m is the apparent magnitude, M is absolute magnitude and di \s the luminosity 
distance which depends upon the geometry of the space and its matter content. In the 
standard scenario a white dwarf accretes mass from a companion star. Once the Chan- 
drasekhar mass limit is reached, the burning of carbon is ignited in the interior of the white 
dwarf. This process propagates to the exterior layers leading to a complete destruction of 
the star. The physics of these objects is not completely understood yet. It requires the 
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Figure 1.1: Lines of constant HqIq in the — plane. From top left to bottom right 
Hoto = (1.08,0.94,0.9,0.85,0.82,0.8,0.67). 

use of numerical simulations from which it appears that the thermonuclear combustion is 
highly turbulent. Such theoretical uncertainties prevent us from having reliable predic- 
tions of possible evolutionary effects. It is also matter of debate as to whether the history 
of the supernova progenitors can have important effects in the final explosion (see [49] for 
a general review) . For these reasons it is a rather unreliable assumption that supernova la 
are perfect standard candles. However the observations show the existence of an empirical 
relation between the absolute peak luminosity and the light curve shapes. There are also 
correlations with the spectral properties. Using such relations it is possible to reduce the 
dispersion on magnitude of each supernova to within 0.17 magnitudes allowing them to 
be used for cosmological distance measurements (see [50] and references therein). The 
Supernova Cosmology Project [51] and the High-Z Supernova Research Team [52] have 
observed and calibrated a large sample of supernovae at low and high redshifts. The re- 
sult of their analysis [52,53] shows that distant supernovae are on the average about 0.20 
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magnitudes fainter than would be expected in a Milne universe (empty). The likelihood 
analysis, due to the deg eneracy of the luminosity distance with the values of ^\ and ^rn: 
constrains these parameters in a region approximated by O.S^lm — O.GJIa ~ —0.2 it 0.1. 
The data give evidence for a non vanishing cosmological constant. Including the farthest 
supernova Sn 1997f with z ~ 1.7 [54], the data analysis shows that at redshift z ^ 1.2 
the Universe was in a decelerating phase. However the presence of possible systematic 
uncertainties has attracted some criticism. In particular there could be a dimming of the 
light coming from the supernovae due to intergalactic dust. Moreover the Sn la might have 
an evolution over the cosmic time, due to changes in characteristics of the progenitors so 
as to make their use as standard candles unreliable. The argument against the extinction 
is that high-redshift supernovae suffer little reddening. While the fact that their spectra 
appear similar to those at low-redshift seems to exclude the possibility of evolutionary 
effects in the data. For instance in [55], the Hubble diagram of distant type la supernovae 
segregated according to the type of host galaxy has been analysed. The results shows 
that host galaxy extinction is unlikely to systematically affect the luminosity of Sn la in 
a manner that would produce a spurious cosmological constant. In reality only a theoret- 
ical prediction, not available at the present time, would convince the entire community. 
Nevertheless none of these systematic errors can reconcile the data with a vanishing A. 

1.4 Cosmic complementarity 

Figure E21 shows the region of the — constrained by different cosmological observa- 
tions. It appears evident that the joint analysis of the recent CMB, large scale structure 
and Sn la observations, previously reviewed, indicate a region of the parameter space 
where they are consistent [56-58] . In particular the likelihood contours of Sn la and CMB 
are orthogonal and therefore their combination breaks the geometric degeneracy between 
Qm and Such consistency tells us that the Universe is nearly flat, the structures we 
observe today are the result of the growth of initial density fluctuations characterized by a 
nearly scale invariant spectrum as predicted by inflationary scenarios. The matter content 
of the Universe consist of baryons (3%), while most of the clustered matter consists of 
cold dark particles that account for only 30% of the total energy density. About 67% of 
the matter is in a 'dark energy' form and is responsible for the present accelerated expan- 
sion. In spite of any astronomical uncertainties, the limits on J7a impose that A < 10~^^ 
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Figure 1.2: 1,2 and 3 a confidence contours on fi^ — plane determined from Sn la, 
CMB and Chandra /gas(z) data (from [45]) 

GeV^. Understanding the origin and the smallness of this term is the challenge of modern 
theoretical physics. In principle nothing prevents the existence of a A-term in the Ein- 
stein equations. However if A appears on the left hand side, the gravitational part of the 
Einstein-Hilbert action will depend on two fundamental constants, G and A, which dif- 
fer widely in scale. For instance the dimensionless combination of fundamental constants 
{Gh/ c?)"^ K < 10^^^^. On the other hand we know that several independent phenomena 
can contribute as an effective A term on the right hand side of the Einstein equations [2] . 
However, as we shall see in the next chapter, in order to reproduce the small observed 
value, these terms have to be fine tuned with bizarre accuracy. Therefore the solution to 
such an enigma may well lead us to the discovery of new physics. 



Chapter 2 

An explanation for the dark 

energy? 



The cosmological constant can be naturally interpreted as the energy contribution of 
the vacuum. However its measured value turns out to be extremely small compared to 
the particle physics expectations. Therefore alternative candidates for the dark energy 
component have been considered. In particular a light scalar field rolling down its self- 
interacting potential can provide the missing energy in the Universe and drive a late time 
phase of accelerated expansion. A lot of effort has gone into justifying the existence of 
this field, called quintessence, within the context of particle theories beyond the Standard 
Model. Different versions of this original idea have been developed in the literature. In this 
Chapter we will review the vacuum energy problem. Then we will discuss the application 
of the anthropic principle to the solution of the 'coincidence problem'. We will describe 
the main characteristic of a minimally coupled quintessence scenario. At the end we will 
analyse the dynamics of two specific scalar field models. 

2.1 Vacuum energy 

It was initially pointed out by Y.B. Zeldovich [59] that in Minkowski space-time, Lorentz 
invariance constrains the energy momentum tensor of zero point vacuum fluctuations to 
be proportional to the Minkowski metric, i.e. T^y'^ = const, x diag{l,—l,—l,—l). This 
relation can be generalized to the case of a curved space-time with metric g^y. The prin- 
ciple of general covariance requires that T^""^ oc g^^, which has the form of a cosmological 
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constant. This implies that in General Relativity, since the gravitational field couples 
through the Einstein equations with all kinds of energy, the vacuum energy contributes 
to the total curvature of space-time. The vacuum state of a collection of quantum fields, 
that describes the known forces and particles, is defined to be the lowest energy density 
state. If we think of the fields as a set of harmonic oscillators the total zero point energy 
is given by: 



that diverges as (ultraviolet divergence). However any quantum field theory is valid up 
to a limiting cut-off scale, beyond which it is necessary to formulate a more fundamental 
description. Consequently the integral Eq. (|2.1j) can be regularized imposing a cut-off 
kmax, and we obtain 



If we set the cut-off kmax at the Planck scale, the energy density of the vacuum is p^ac ^ 
(10^^ GeV)^ which is about 120 orders of magnitude larger than the observed value of p\. 
Fixing the cut-off scale at the QCD phase transition, k^ax = ^qcd, we find pvac^ ~ 10""^ 
GeV^ which is still 44 orders of magnitude above the expected one. On the other hand 
if Supersymmetry is realized in nature, the cosmological constant vanishes because the 
vacuum energy contribution of the bosonic degree of freedom exactly cancels that of the 
fermionic ones. However, because we do not observe super-particles, Supersymmetry must 
be broken at low energy. This implies that the cosmological constant vanishes in the early 
Universe and reappears later after SUSY breaking. Assuming that Supersymmetry is 
broken at MsusY ~ ITeV the resulting p^ac is about 60 orders of magnitude larger then 
the observational upper bounds. Hence any cancellation mechanism will require a bizarre 
fine tuning in order to explain the huge discrepancy between pyac and pA- By the present 
time we do not have any theoretical explanation for this cosmological constant problem. 
Moreover such a tiny value presents an other intriguing aspect. In fact we could ask why 
A has been fixed at very early time with such an extraordinary accuracy that today it 
becomes the dominant component of the Universe. In other words we should explain why 
the time when A starts dominating nearly coincides with the epoch of galaxy formation. 
This is the so called coincidence or 'why now' problem. The solution to the cosmological 
constant problem will provide an explanation also for this cosmic coincidence. On the 
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other hand it could be easier to justify a vanishing cosmological constant assuming the 
existence of some unknown symmetry coming from quantum gravity or string theory [60] . 
As we shall see in the following sections, alternative scenarios of dark energy formulate 
the initial condition problem and the coincidence problem in a different way. 

2.2 Anthropic solutions 

The use of anthropic arguments in cosmology has been often seen as an anti-scientific 
approach. However a different use of the 'Anthropic Principle' has been recently proposed 
in the literature and for a review of the subject we refer to [61]. We should always have 
in mind that at the speculative level our Universe can be one particular realization of 
possible universes. Therefore the fact that we live in this Universe makes us privileged 
observers, since under other circumstances we would not be here. For instance several 
authors pointed out that not all values of A are consistent with the existence of conscious 
observers [62-64]. The reason is that in a flat space-time the gravitational collapse of 
structure stops at the time t tjy, as consequence universes with large values of A will not 
have galaxies formed at all. This argument can be used to put an anthropic bound on pA 
by requiring that it does not dominate before the redshift z^ax when the earliest galaxy 
formed. In [64] assuming z^ax = 4 it was found p\ < p^. However it was suggested 
in [65, 66] that observers are in galaxies and therefore there is a conditional probability 
to observe a given value of A. In particular this value will be the one that maximizes the 
number of galaxies. In such a case the probability distribution can be written as 



where ^*{pa) is the a priori probability density distribution and ^{ph) is the average 
number of galaxies that form per unit volume with a given value of pA- The calculation 
of v{pa) can be done using the Press-Schechter formalism. Assuming a flat a priori 
probability density distribution the authors of [67] found that the peak of 3^{ph) is close 
to the observed value of A. This anthropic solution to the cosmological constant problem 
would be incomplete without an underlying theory that allows A to take different values 
and predicts a flat ^*{pa)- The recent developments in string/M theory seems to provide 
a natural framework where such issues can be addressed (see [68]). 
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2.3 Quintessence 

A non-anthropic solution to the cosmic coincidence problem would be an exotic form of 
matter playing the role of dark energy. The existence of such a component should be the 
prediction of some fundamental theory of particle physics. For instance it was initially 
suggested that a network of topological defects could provide such a form of energy [69-71]. 
In fact topological defects are characterized by a negative value of the equation of state 
wx ^ — 1/3 and lead to an accelerated expansion if they dominate the energy budget of 
the Universe. However these models are ruled out by current cosmological observations. 
On the other hand, long before the time of Sn la measurements, it was considered that an 
evolving scalar field, called quintessence, could take into account for the missing energy 
of the Universe [72-78]. In this scenario the cosmic coincidence problem is formulated 
in a different way. In fact the evolution of the quintessence is determined by the initial 
conditions and by the scalar field potential. Consequently there would be no coincidence 
problem only if the quintessence becomes the dominant component today independently of 
the initial conditions, that have been set at very early time. It was pointed out by Zlatev, 
Wang and Steinhardt [79, 80] that viable quintessence potentials are those manifesting 
'tracking' properties. In these cases, for a wide range of initial conditions, the scalar field 
evolves towards an attractor solution such that at late time it dominates over the other 
matter components. However such a time will depend on the energy scale of the potential 
and is fixed in way such that pQ reproduces the observed amount of dark energy. In other 
words the tracker quintessence solves the initial conditions problem, but the 'why now' 
problem is related to the energy scale of the model. If such a scale is consistent with 
the high energy physics scales there is no fine tunning and the fact that the acceleration 
starts only by the present time does not have any particular meaning. On the contrary 
the coincidence problem would result if such a scale is much smaller than any particle 
physics scale, because this will require a fine tuning similar to the cosmological constant 
case. As we shall see consistent quintessence model building is a difficult challenge [81]. 
Cosmic coincidence is absent in quintessence models where the scalar field is non-minimally 
coupled to the cold dark matter [82-85]. Coupling with baryons is strongly constrained 
by tests of the equivalence principle, however a coupling with cold dark matter cannot be 
excluded. In this case the coupling will naturally produce the gravitational collapsing time 
scale of the order of the time when the Q field starts dominating, to ^ tq. Moreover in 
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these models structure formation can occur even during the accelerated phase of expansion 
and consequently no coincidence have to be explained [86]. It may be argued that for this 
class of models to fully succeed what has to be explained is the strength of the couplings. 
However the non universality of the couplings may arise in the context of brane models, 
where dark energy and dark matter belong to an hidden sector. 

2.3.1 Scalar field dynamics 

A multiple fluid system consisting of a scalar field, pressureless matter and radiation 
interacting through the gravitational field is described by the action: 

S = J d^x^gR + j d^x^i^Q + ^m), (2.4) 

where R is the Ricci scalar and Lagrangian density of the matter and radiation and 
is the Lagrangian density of the quintessence field which is given by: 

^Q = ^d^^Qd^Q-V{Q). (2.5) 

The scalar field energy-momentum tensor then reads as 

Tq^^ = d^Qd.Q - g^, Qs^ga^Q - V{Q)^ . (2.6) 

In a FRW flat Universe for a nearly homogeneous scalar field, the quintessence pressure 
and energy density are pg = — V and pq = + V. The quintessence behaves as 
perfect fluid with a time dependent equation of state which is given by: 

^^^ft^XM. (2.7) 

The scalar field evolution is described by the Klein-Gordon equation, 

Q + 3ifQ + ^ = 0, (2.8) 
dQ 



with 

8ttG 



n2 

Pm + Pr + ^ + V{Q) 



(2.9) 



3 

where pm and pr are the matter and radiation energy densities and evolve according to 
the energy conservation equation Eq. ()1.7() . As an example we analyse the dynamics of 
this system in the case of an Inverse Power Law potential [73,79]: 

Aa+4 

V{Q) = (2.10) 
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where a = 6 and we set A such that Qq = 0.7. We solve numerically the equations of 
motion. In figure lTTI we show the evolution with redshift of the energy density of radiation 
(blue dash line), matter (green dot line) and quintessence (red solid line). Different red 
lines correspond to different initial conditions, we may distinguish two distinct behaviours. 
For initial values of the quintessence energy density larger than matter energy density, 
/3q > pq rapidly decreases. This regime is called kination. In fact, as we can see in 
figure HTIh . the kinetic energy rapidly falls off while the potential energy remains nearly 
constant. The overshooting is then followed by the frozen field phase, where the energy 
density is dominated by the potential. During this period pq remains constant until the 
kinetic energy becomes comparable with the potential one and the field reaches the tracker 
regime. In the tracker solution the kinetic and potential energies scale with a constant 
ratio, therefore the equation of state is constant or slowly varying. For a given potential the 
existence of a period of tracking is guaranteed by the condition F = V"V/{V')'^ > 1 [80]. 
Moreover during this phase the value of the quintessence equation of state mimics the 
value wb the background component according to the relation: 



For inverse power law potentials F = 1 + a~ . At late time the field leaves the tracker 
solution, the potential energy starts dominating over the kinetic one (figure [T^)) and the 
equation of state tends to negative values. During this final phase the field becomes the 
dominant component of the Universe and drives the accelerated expansion. The same 



quintessence starts its evolution in the frozen regime. Such a behavior occurs over a 
range of initial conditions that covers more then 100 orders of magnitude, consequently 
the quintessence dominated period is obtained with no need of fine tunning of the initial 
conditions. For tracker models the final value of the equation of state depends on the 
parameter A or equivalently on fig and on the slope of the potential. In general for large 
values of Qq or fiat potential wq —1, in the case of the Inverse Power Law potential 
small values of a corresponds to large negative value of Wq (figure IT^ . 

2.3.2 Quintessential problems 

The existence of the tracker phase cancels any knowledge of the initial condition of the 
scalar field providing an elegant way of solving the coincidence problem. Nevertheless it is 




WB - 2(F - 1) 




arguments hold for initial conditions corresponding to pg < p^. In such a case the 
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Figure 2.1: Energy density versus redshift for radiation (blue dash line), matter (green 
dot line) and quintessence (red solid line). As we may notice for a large range of initial 
conditions the quintessence energy density converges to the tracker solution. 



a rather difficult task to build consistent particle physics models of quintessence. One of 
the reasons it that for a given tracker potential, V{Q) = Af(Q), the cosmic coincidence is 
fully solved only if the scale A is consistent with the energy scale of the underlying particle 
physics theory that predicts the shape f{Q) of the quintessence potential. For example let 
us consider the Inverse Power Law model. It was shown in [87, 88] that it can be derived 
from a Supersymmetric extension of QCD. By the present time the condition |V^V^I < 1 
has to hold in order to guarantee the Universe is accelerating, this implies that today the 
scalar field Q ~ Mpi. Since the observations suggest pg ~ pc, for values of the slope a > 6 
we find that A ~ 4.8 x 10^ GeV, a very reasonable scale from the particle theory point of 
view. On the other hand for tracker models the slope of the potential is constrained by 
measurements of the present value of the quintessence equation of state Wq that indicate 
a low value of a. But, as we can see in figure 12.41 for values of a < 6 the energy scale 
is much smaller than any known particle physics scale. Therefore the Inverse Power Law 
seems not to be a viable quintessence model. Alternative models have been proposed in 
the literature, they can be distinguished into two categories, dilatonic and supersymmetric 
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Figure 2.2: Evolution of the kinetic (blue dash line) and potential energy density (red 
solid line) at early times (a) and after matter-radiation equality (b). 
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Figure 2.3: Equation of state versus redshift for an Inverse Power Law potential with 
a = 6 (blue solid line) a = 4 (red solid line) and q = 2 (green solid line). 
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Figure 2.4: Energy scale versus a for the quintessence Inverse Power Law potential. 

quintessence. The former class of models assume that the quintessence field is the dilaton, 
this possibility has been studied in [89,90]. The dilaton is predicted by all string theory 
models and it couples to all the fields including gravity [91]. Therefore it could be a 
good candidate for dark energy. However it predicts a running of the different coupling 
constants, that are strongly constrained by present observations and the violation of the 
equivalence principle. Nevertheless these models deserve more investigation and some of 
these issues have been recently addressed in [92,93]. It is worth mentioning that the 
non-minimally coupled scalar field models which solve the coincidence problem belong to 
this category [82-85,94,95]. In the second class of models the quintessence is one of the 
scalar fields predicted by Supersymmetric extensions of the Standard Model of particle 
physics. In particular a lot of effort has been recently devoted to the formulation of 
viable quintessence models in the context of Supergravity theory. In fact it was noticed 
that Inverse Power Law potentials generated by Supersymmetric gauge theories are stable 
against quantum and curvature corrections, but not against Kahlerian corrections [96-98]. 
Since Qo ~ Mpi, Supergravity (SUGRA) corrections cannot be neglected and therefore 
any realistic model of quintessence must be based on SUGRA. It is argued that also this 
class of models lead to violation of the equivalence principle, the quintessence field can in 
fact mediate a long range fifth force that we do not observe. However we want to stress 
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that in the context of Supersymmetric theories this is not the main problem. The reason 
is that the quintessence can belong to a hidden sector of the theory that couples only 
gravitationally to the visible sector. It was pointed out in [81] that any supersymmetric 
inspired model has to address two specific issues. The first one concerns the case of 
Supersymmetry breaking quintessence, if the quintessence field belongs to a sector of the 
theory that breaks Supersymmetry, because of the shape of the potential it turns out 
it cannot be the main source of breaking. In such a case SUSY can be broken by the 
presence of an F-term that leads to an intolerably large vacuum energy contribution that 
completely spoils the nice properties of the quintessence potential. The other difficulty 
arises from the coupling of quintessence to the field responsible for the supersymmetry 
breaking. Such a coupling leads to corrections of the scalar field potential such that the 
quintessence acquires a large mass. Some alternatives have been recently investigated, for 
instance a way of avoiding such problems has been considered in [99], where a Goldstone- 
type quintessence model in heterotic M-theory has been proposed. 



2.4 Supergravity inspired models 

We now review some of the properties of two models proposed in the context of Supergrav- 
ity theories: the Exponential Times Inverse Power Law potential and the Two Exponential 
potential. 

2.4.1 Exponential Times Inverse Power Law potential 

The authors of [96] have shown that taking into account Supergravity corrections to the 
Inverse Power Law potential, the quintessence potential takes the form: 

V{Q) = ^^^^ > (2-12) 
where k = l/Mp^. This potential is an improvement the Inverse Power Law. In fact the 
dynamic remains unchanged during the radiation and the matter dominated era, while 
the presence of the exponential term flatten the shape of the potential in the region 
corresponding to the late time evolution of the scalar field. This allows for more negative 
values of the equation of state today independently of the slope of the inverse power law. 
Consequently we can have a reasonable particle physics energy scale even for large values 
of a. For instance for a = 11 and Qq = 0.7 we have A ~ 10^^ GeV and the present value 
of the equation of state is Wq = —0.82, in better agreement with observational constraints. 
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2.4.2 Two Exponential potential 

The dynamics of cosmologically relevant scalar fields with a single exponential potential 
has been largely studied in the literature and within a variety of contexts (see for instance 
[77]). The existence of scalar field dominated attractor solutions is well known, however 
for this class of models an accelerated phase of expansion can be obtained with an extreme 
fine tunning of the initial conditions [102]. It has been shown in [100,101] that quintessence 
Supergravity inspired models predict the scalar field potential to be of the form: 

V[Q) = Ml.1 (^e"v^(^-^) + e^v^(^-^)) , (2.13) 

where A is a free parameter, while B has to be fixed such that Mp^e~^^ ~ pg. This 
potential has a number of interesting features. As it has been pointed out by the authors 
of [100] in the form given by the Eq. ()2.13|) all the parameters are of the order of the 
Planck scale. Only B has to be adjusted so that Mpie~^^ ~ pg, it turns out to be 
B = ^(lOO)Mp;. For a large range of initial conditions the quintessence field reaches the 
tracker regime during which it exactly mimics the evolution of the barotropic fluid and at 
some recent epoch it evolves into a quintessence dominated regime. It is useful to rewrite 
the two exponential potential as: 

V{Q) = (^e-"«/^^^' + e-^^/^^^') , (2.14) 

where M is the usual energy scale parameter. For a given value of 0,q the slopes a, [3 fix 
the final value of the equation of state. The sign of the slopes distinguish this class of 
potentials into two categories: a > /3 > and a > 0, /3 < 0. For both the cases the Q-field 
initially assumes negative values and rolls down the region of the potential dominated by 
the exponential of a. When it reaches the tracker regime its equation of state exactly 
reproduces the value of the background dominant component, wq = wb and the energy 
density is given by 

= 3{WB + \)lo?. (2.15) 

The late time evolution is determined by the value of /? that fixes the present value of 
the equation of state. In the case of slopes with the same sign w^q — > — 1 for small value 
of /3, while for a and /3 with opposite sign the scalar field reaches by the present time 
the minimum of potential at Qmin/Mpi = ln{—a/ [5) / {a — /3). Consequently uIq — 1 
after a series of small damped oscillations. This can be seen in figure l231 where we plot 
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Figure 2.5: Evolution of the quintessence energy density and equation of state for param- 
eters {a, (3): blue solide line (20,0.5); red dashed line (-20,-20) and Qq = 0.7. 

the evolution of the quintessence energy density parameter 0,q and the equation of state 
wq for a = 4, /3 = 0.02 and a = 20, f3 = —20. In the latter model the equation of 
rapidly drops to —1 after few damped oscillations, while the former shows a more smooth 
behaviour. It is worth remarking that for the two exponential potential with same sign of 
the slopes the accelerated phase can be a transient regime. This can occurs for large value 
of (3, in such a case after a short period dominated by the potential energy the scalar field 
acquires kinetic energy so that the equation of state can be wq > —1/3. As we may note 
in figure because Eq. (|2.15j) holds during the tracker regime, for small values of a the 
quintessence energy density assumes non-negligible values at early times. Such an early 
contribution during the radiation dominated era is constrained by nucleosynthesis bound 
to be ^q{1 MeV) < 0.13. This implies that a > 5.5. Such a limit has pushed toward 
larger value by a new analysis of the Big-Bang nucleosynthesis [103] that constrains ^q{1 
MeV) < 0.045 at 2a. In principle this bound can be avoided if the tracker regime starts 
after Big-Bang nucleosynthesis. However this can be obtained only by tuning the scalar 
field initial conditions in a restricted range of values. Moreover a recent analysis of the 
large scale structure data and CMB measurements strongly constrain the value of 
during the matter dominated era [104]. 




Chapter 3 

Quintessence field fluctuations 



The cosmological constant is a smooth component and therefore does not play any active 
role during the period of structure formation. On the contrary a peculiar feature of the 
quintessence field is that it is spatially inhomogeneous just as any other scalar field. There- 
fore it is possible that the clustering properties of the dark energy can play a determinant 
role in revealing the nature of this exotic fluid. In this Chapter we introduce the scalar 
fleld fluctuation equations in the Newtonian gauge. We derive an analytical solution of 
the quintessence perturbation in the tracking regime and describe the behaviour during 
different cosmological epochs. We then present the numerical analysis of the perturbations 
in a multiple fluids system in the particular cases of an Inverse Power law potential and 
the two exponential potential. 



3.1 Quintessence perturbations in Newtonian gauge 

The evolution of minimally coupled scalar fleld perturbations has been studied in a number 
of papers. For instance in [73, 77, 105, 106] the analysis has been done in the synchronous 
gauge, while in [107-109] the authors have used the Newtonian gauge. In what follows 
we use the Newtonian gauge. When compared to the synchronous one it has a number of 
advantages. In fact since the gauge freedom is fully flxed there are no gauge modes that 
can lead to misleading conclusions about the evolution of superhorizon modes. Besides, 
the metric perturbations play the role of the gravitational potential in the Newtonian 
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limit. The equations we need to linearize are the Einstein's equations 

Rtiu - ^9fiiyR = SvrGT^i., (3.1) 
the Klein-Gordon equation 

^ d,{V^g^''dM) + ^ = o (3.2) 



and the conservation equation of the stress energy tensor of the different matter compo- 
nents 

= 0, (3.3) 

In the Newtonian gauge the line element in a spatially flat FRW background reads as 

ds^ = (1 + 2$)dt2 _ a2(t)(l - 2<^)dx'dxi, (3.4) 

where $ is the metric perturbation and t is the real time. We consider a multiple fluid 
system composed of a scalar field, cold dark matter and radiation. Expanding the fluid 
variables at first order around the homogeneous value we have: 

pl = Pi{l + 6i), (3.5) 

where Si is the density perturbation of the i-th component, 

Q^ = Q + SQ, (3.6) 

where 5Q is the scalar field fluctuation and Q is the homogenous part of the quintessence 
field. The linearized Eq. (|3.H) . Eq. (|3.2j) and Eq. (|3.3jl provide a set of differential equations 
for the metric, scalar field, radiation and cold dark matter perturbations. In Fourier space 
the equations are: 

5Q + 3HSQ + ^6Q + —rSQ = 4Q^-2—<^, (3.8) 

aQ^ dQ 

5c = -Vc + 3$, (3.9) 
a 

K = (3.10) 
a 

Ah 

8r = —Vr + A^, (3.11) 

6a 



3.2 Evolution of perturbations 



28 



Vr = —^6r - (3.12) 

4a a 

where 5c and 6r are the density perturbations of cold dark matter and radiation, while Vc 
and Vr are the corresponding velocity perturbations. The perturbed quintessence energy 
density and pressure are, 

5pq=pq5q = Q6Q-Q^<!> + ^6Q, (3.13) 
5pQ = Q6Q - _ ^§Q, (3.14) 

3.2 Evolution of perturbations 

3.2.1 Analytical solution in the tracker regime 

Prom Eq. (|3.8|) we note that the perturbations of the other fluids feed back onto the scalar 
field perturbations through the gravitational potential Deep during the radiation and 
the matter dominated eras the gravitational potential is constant and therefore the first 
term on the right-hand-side can be neglected. As a first approximation we ignore the 
second term as well, hence Eq. (|3.8|) becomes 

fV' 

We can find an analytical solution during the tracker regime when 

%~A^,H^ ,3.16) 

with Aq a constant that depends on the specifics of the potential and on the equation of 
state of the background dominant component. This can be obtained as follows. Consider 
the adiabatic definition of the sound speed associated with the quintessence field: 

= 1 + ^^. (3.17) 

mq 

During the tracker regime the quintessence equation of state is nearly constant implying 
that Cq « WQ = const. Hence by differentiating Eq. (|3.17j) with respect to time and 
dividing by H we have 

'l + ;^-^(34 + 5) , (3.18) 



5Q + m5Q +[ — + -T7^]5Q = ^- (3.15) 



c: 
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From Eq. (|3.18|) using the second of the Hubble equations H / = —3(1 + wb)/'^, with 
wb being the equation of state of the background dominant fluid, we finally obtain: 

^ = \{l-cl){wB + cl + 2)^Al. (3.19) 

It is useful to rewrite Eq. (|3.15|) in conformal time, which is defined as d/dr = a{t)d/dt, 

6Q" + 2J^6Q' + (fc2 + a^H^A^) 6Q = 0, (3.20) 

where a prime denotes the derivative with respect to r. In the radiation dominated era 
the scale factor evolves as a = t/2, while the Hubble rate is given by H = 2/r, hence 
Eq. (|3.2flj) becomes: 

5Q" + 26Q' + (k^ + AA'q)5Q = 0, (3.21) 

with Aq being the value of Aq in radiation dominated era. During the tracker regime 
Cq ^ Wr = 1/3 and A^q ~ 4. In such a case the characteristic roots of Eq. 1)3.21(1 are 
complex and the solutions are of the form: 

6Q = e~^^ (Ci cos VrT + C2 sin z/^r) , (3.22) 

where Ci and C2 are integration constants and 



fc2 + 4A^-l. (3.23) 

Eq. (|3.22() is a damped oscillatory solution with frequency m^. On large scales [k < 1 
Mpc~^) the frequency of these oscillations is scale independent and is set by the tracker 
properties (^g)- On the other hand if A^q = 0, Eq. (|3.21() has real roots and the solutions 
contain a constant mode and a decaying one. Similar solutions can be found in the matter 
dominated era. The addition of the source term $ on the right-hand-side of Eq. (|3.<S|) 
leads to an attractor solution for the quintessence perturbations in the long-wavelength 
limit, 6Q{t) SQc- In fact from Eq. (|3.8() we obtain: 

6Qc ~ -2^<^c, (3.24) 

that is constant as long as the quintessence is in the tracker regime. As a consequence 
this solution does not hold in the kinetic and particularly in the potential phase when the 
quintessence exits away from the tracker. 
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Figure 3.1: Evolution of the scalar field perturbation 6Qk for k = 0.00001 (blue solid line), 
k = 0.001 (red dash line) and A; = 0.1 (green dash line). 

3.2.2 Numerical analysis 

In this section we present the analysis of the numerical solution of the system of equa- 
tions Eq. H3.7I3.12() . where we have imposed adiabatic initial conditions. We consider a 
quintessence SUGRA inspired model described in Chapter El The scalar field potential 
is specified by Eq. (|2.12|) where we set the slope a = 6 and the parameter M such that 
today = 0.7. The initial conditions for the homogeneous part of the scalar field Q 
have been set such that the tracker solution is reached deep in the radiation dominated 
era. In figure 13.11 we plot the behaviour of 5Q against the conformal time for three dif- 
ferent wavenumbers. We can see that as the system enters the tracker regime, decaying 
oscillations are set for all the three modes with a constant frequency given by Eq. ()3.23|1 . 
Since A;^ is negligible in Eq. (|3.23j) the frequency of these oscillations is the same for all 
modes. When the second term in the right-hand-side of Eq. ()3.8|) becomes comparable 
to the term V^qq the fluctuations evolves onto the attractor solution given by Eq. ()3.24p . 
In the specific case we consider there is not an exact tracking and consequently the ratio 
V^q/V^qq scales linearly and not steadily as in the exact tracking case. This explains why 
5Q increases with time in the long wavelength mode. As different scales cross the horizon 
they leave this attractor solution. For instance the shortest wavelengths {k = 0.1) enter 
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Figure 3.2: Phase diagram for long wavelengtli quintessence fluctuations with two different 
initial conditions. The attractor point 5Qc (Eq. I3.24|) is a transient attractor. 

the horizon after matter radiation equality. Due to the decay of the gravitational poten- 
tial the first term on the right-hand-side of Eq. (|3.8j) becomes large leading to a decaying 
solution with small damped oscillations whose frequency is set by the effective mass term 
k on the left-hand-side of Eq. 1)3. 8|) . The presence of an attractor solution for the scalar 
field perturbations is evident in figure where we plot the phase space 6Q — 6Q. As the 
system evolves in the tracker regime, 6Q and 6Q spiral toward the attractor point dQc- 
As soon as the tracker terminates the long wavelength fiuctuations have converged to the 
same values and their evolution is indistinguishable. In minimally coupled quintessence 
models the effects of the scalar field perturbations on the structure formation process are 
negligible. This is because the quintessence interacts only gravitationally with the other 
fluid components and it is usually a subdominant component during the epochs that are 
relevant for the formation of the early structures. Therefore the clustering of the dark 
energy can have detectable effects only at late times through the time evolution of the 
gravitational potentials. As we shall see in Chapter IHl this produces a characteristic im- 
print in the CMB anisotropy power spectrum. In figure 13.31 we plot the behaviour of the 
gravitational potential $ normalized to the cosmological constant case for k = 0.0004 (up- 
per panel) and k = 0.04 for three different quintessence models whose dynamics have been 
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described in Cliapter |2 (model A) the SUGRA potential Eq. (|2.12|) with slope a = 6 
(blue solid line); (model B) the two exponential potential Eq. (|2.1I'{|) with a = 20 and 
(3 = —20 (red dash hne) and (model C) a = 4 and /3 = 0.02 (blue dash dot line). We 
may note a number of features that arises from the presence of dark energy perturbations 
and a different expansion rate of the Universe compared to the ACDM model. Let us 
consider the mode k = 0.0004 that crosses the horizon during the matter era. For model 
B (red dash line) the quintessence energy density is always negligible during the matter 
era and consequently there are no effects due to the perturbations in the evolution of 
the gravitational potential. On the other hand because of the scalar field dynamics the 
late time expansion rate of the Universe is different compared to the A case, hence <^ /<I>a 
diverges from the unity. In contrast in models A and C the quintessence energy density 
dominates earlier, especially in model C where Qq is a large fraction of the total energy 
density (see figure IT^ . In this case we can distinguish two different evolutionary regimes 
of the gravitational potential. At redshift z > 10 the Universe is matter dominated and the 
quintessence is in the tracker solution with an equation of state very close to the matter 
value. Therefore during this period the expansion rate is the same as in the A model and 
we would expect ^/^a ~ 1- The fact this ratio deviates from unity is due to the presence 
of the dark energy perturbations. On the other hand the evolution of the gravitational 
potential at z << 10 is caused by the accelerated phase of expansion. This can be seen 
by the change of the slope of ^/^a at a redshift 2; ~ 2 when the acceleration starts. The 
late time decay of the gravitational potential sources the formation of CMB anisotropy 
through the late Integrated Sachs- Wolfe effect (ISW). The same arguments hold for the 
mode k = 0.04 that enters the horizon soon after matter-radiation equality. However it is 
worth noticing that in the case of model C, the early contribution of dark energy causes 
the scalar field perturbations to produce a bigger effect on the decay of leading to a 
larger early ISW effect. From this analysis we conclude that different dark energy models 
will lead to a characteristic signature in the CMB power spectrum and will prove this 
point in a more general way in Chapter El 
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Figure 3.3: Evolution of tiie gravitational potential normalized to the cosmological 
constant case for a SUGRA model (blue solid line), two exponential potentials with slopes 
(—20,20) (red dash line) and (4,0.02) (green dash line) respectively. The top panel shows 
the mode with k = 0.0004 and the lower panel shows k = 0.04. 



Chapter 4 

Constraining the quintessence 

potential 



Cosmological distance measurements test the expansion rate of the Universe at different 
redshifts. Therefore they can be used to constrain the dark energy properties and even- 
tuahy disprove the ACDM model. A lot of effort has gone into constraining the value 
of the dark energy equation of state with the current available data. In this Chapter we 
will review some of the constraints obtained in previous works. Then we will describe the 
analysis of the Sn la data and the measured position of the CMB peaks that we performed 
to constrain the shape of a general quintessence potential [110]. We will comment on our 
results and compare them to those obtained in other related work. 



4.1 Upper bounds on the cosmic equation of state 

Different dark energy models are usually constrained assuming that the dark energy be- 
haves as a perfect fluid with a constant equation of state wx = Px/px ^ —1- This can 
take into account for models such as a network of topological defects with wx = —n/3 (n 
is the dimension of the defect) [111] and as first approximation time dependent tracker 
quintessence models with wq > — 1. Different values of wx will lead to a different expan- 
sion rate and consequently to a different cosmological distance. In particular for larger 
negative values of wx the phase of acceleration starts at earlier times and therefore cosmo- 
logical distance indicators will appear farther than in models with smaller negative values 
of Wx- This can be seen in figure UHl we plot in the solid line the redshifts Zacc when the 
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Figure 4.1: Redshifts Zacc (solid line) and Zeq (starred line) corresponding to the start 
of the accelerated phase and the dark energy dominated epoch versus a constant wx for 
VLx = 0.6,0.7,0.8 (in blue, red and green lines respectively). 

Universe starts accelerating and in a starred line the redshift Zde when the dark energy 
starts dominating as a function of wx for different values of Ox- For a flat Universe The 
redshift Zacc is given by 

(1 + Zaccf = -(1 + ^wx)-r^^, (4.1) 



and the redshift z^^, is 







-l/wx 



Note for wx < —0.6 we have Zacc > z^e- For a flat Universe the effect on cosmological 
distances produced by varying wx is degenerate with 0^ = 1 — Vtx and with the Hubble 
parameter Hq. By making use of the Sn la data the authors of [53, 112] find the constraint 
Wx < —0.55 (2(t) for > 0.1. We will comment in the next Chapter about the bias effect 
introduced in the analysis of data by assuming a constant equation of state. The position 
of the Doppler peaks in the CMB anisotropy power spectrum is another indicator of the 
cosmological distance to the last scattering surface and hence can be used to constrain 
wx- In [113] a constant effective equation of state 
J wxia)^xia)da 
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Figure 4.2: Upper panel: relation between Weff and Om. for an Inverse Power Law po- 
tential with slope p and an Exponential potential V oc exp(mp//Q). Lower panels: One 
dimensional likelihoods as a function of 0,^ for CMB+LSS (dot) and Sn la (solid) data 
(from [113]). 

has been constrained with the Sn la, CMB and the LSS data. The authors find Weff < 
—0.65 (2(t) with Weff = —1 being best fit. They comment that tracker models can 
marginally accommodate the Sn la and LSS constraints at the same time (see figure . 
Such an inconsistency can be explained as follows. The supernova luminosity distance 
is sensitive to the value of Weff, while large scale structure data are more sensitive to 
the amount of clustered matter 0^ through the matter power spectrum. This implies 
that dark energy models with different values of w^jf will not change the likelihood in 
Qm of the large scale structure data. In fact assuming the value of the effective equation 
of state to be w^jf > —1, this shifts the Sn la likelihood towards smaller values of 
{i.e. larger values of i^x)- As we will show later, extending the likelihood analysis to the 
slope of the tracker potential enables us to find models that simultaneously fit the whole 
data. Similar constraints from Sn la and CMB data have been found by Efstathiou [114], 
namely wx ^ —0.6 {2a). In a different approach Saini et al. [115] have parameterized 
the luminosity distance and by constraining their parameter space with the Sn la data 
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they have been able to reconstruct the redshift dependence of the dark energy equation of 
state. They find a time dependence with — 1 < wx < —0.86 today and — 1 < wx < —0.6 
at z = 0.83. Waga and Frieman have constrained the slope of the Inverse Power Law 
potential a (see Eq. ()2.10|) ) by making use of the Sn la and lensing statistic data [116]. 
Imposing > 0.3 they obtain an upper limit on the present value of the equation of 
state Wq < —0.67 that implies a < 1.8. The combined analysis of the CMB power 
spectrum measured by the Boomerang and MAXIMA experiments with Sn la and large 
scale structure, limits a constant equation of state to be wx < —0.7 with wx = —1 being 
the favoured value [117]. 

4.2 Parameterized quintessence potential 

We have seen in Chapter [21 that for tracker models the present value of the equation 
of state depends of the shape of the scalar field potential. In particular for SUGRA 
inspired models [96, 100, 101, 118] the equation of state parameter varies in the range 
— 1 ^ Wq < —0.8, while Inverse Power Law potentials require larger values. A general 
potential which can accommodate for these classes of scenarios is: 

^Q) = ^ei('^«)^ (4.4) 

where k = ^/SttG and M is fixed in such a way that today pq = Pc^q, where pc is the 
critical energy density. For (3 = Ea. H4.4|) becomes an inverse power law, while for (3 = 2 
we have the SUGRA potential proposed by [96]. For a = 0, /? = 1 and starting with a large 
value of Q, the Quintessence field evolves in a pure exponential potential [77]. We do not 
consider this case further since it is possible to have a dark energy dominated universe, 
but at the expense of fine tuning for the initial conditions of the scalar field. Larger 
values of (3 mimic the late time evolution of the model studied in [101]. For a,(3^0 the 
potential has a minimum, the dynamics can be summarized as follows: for small values 
of (3 and for a large range of initial conditions, the field does not reach the minimum by 
the present time and hence > —1. For example, if the Quintessence energy density 
initially dominates over the radiation, the Q field quickly rolls down the inverse power law 
part of the potential eventually resting in the minimum with wq ~ — 1 after a series of 
damped oscillations [119]. This behaviour however requires fine tuning the initial value of 
Q to be small. On the other hand, this can be avoided if we consider large values of a 
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Figure 4.3: In (a) the evolution of wq against the red-shift is plotted for different values 
of a and /3. In (b) the behaviour of the deceleration parameter, q, is plotted against the 
red-shift. The acceleration starts {q < 0) earlier for models with an equation of state close 
to that of a true cosmological constant. 

and P (figure Unt)- In these models the fractional energy density of the quintessence field, 
rig, is always negligible during both radiation and matter dominated eras. In fact, for 
small initial values of Q, V{Q) acts like an inverse power law potential, hence as Q enters 
the scaling regime its energy density is sub dominant compared to that of the background 
component. Therefore nucleosynthesis constraints [103] are always satisfied and there are 
no physical effects on the evolution of the density perturbations. The main consequence 
is that for a different value of Wq the Universe starts to accelerate at a different redshift 
(figure 14.3b ') . This implies that different values of a and /? lead to a different luminosity 
distance and angular diameter distance. Consequently by making use of the observed 
distances we may in principle determine an upper limit on Wq, potentially constraining 
the allowed shape of the quintessence potential [120] . 
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4.3 CMB peaks 

The CMB power spectrum provides information on combinations of fundamental cosmo- 
logical parameters. The physical processes responsible for the anisotropy are well under- 
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stood allowing us to accurately predict the shape of the anisotropy power spectrum for 
a given cosmological model. The presence of acoustic oscillations of the photon-baryon 
plasma before the recombination epoch establishes a multiple peaks pattern in the power 
spectrum. For a review of the subject we refer to [121]. During the radiation dominated 
era the equation describing the effective temperature fluctuation AT of the CMB is of the 
form [122]: 

Ar" + c^Ar = 2^'", (4.5) 

where ^ is the gravitational potential sourced by the energy density perturbations in the 
different matter components, Cg = 1/3(1 + "ipb/ Pr) is the sound speed of the barotropic 
fluid and the primes are derivatives with respect to /cr, r being the conformal time and 
k the comoving wavenumber. Eq. ()4.5|) has the form of a perturbed harmonic oscillator, 
therefore when the photons decouple from the baryons their energy carries an imprint 
of such oscillations. The characteristic frequency of these oscillations is fixed by the 
size of the sound horizon at the decoupling = j Cgdr. Therefore we have a series 
of compressions and rarefactions at scales km^sh = "i'^- Today such scales appear at 
angles that are multiple integers of the angular size of the sound horizon at the decoupling 
Gsh = i^sh/ Dk{tis), where Di^{tis) is the distance to the last scattering surface for a space- 
time with curvature k. As a consequence of these the position of the 'Doppler' peaks in 
the power spectrum depends on the geometry of the Universe. For a flat Universe the 
peaks will appear at the multipoles 

Im = mlsh = ( — - 1 ) , (4.6) 

where Cg is the mean sound speed and tq, Tig are the conformal time today and at last 
scattering respectively. However the acoustic oscillations are perturbed by the evolution of 
the gravitational potential ^ on the right-hand-side of Eq. 1)4. 5() , which shifts the position 
of the peaks by an amount that depends on the cosmological parameters that are relevant 
before the recombination. This results in a better estimate for the peak positions being 
given by: 

Im = hh{m - 61 - 5lm), (4.7) 

where 51 is an overall shift [122] and 6lm is the shift of the m-th peak. These corrections 
depend on the amount of baryons Q.}jh?^ on the fractional quintessence energy density at 
last scattering (JIq) and today (^^q), as well as on the scalar spectral index n. Analytic 
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formulae, valid over a large range of the cosmological parameters, have been provided to 
good accuracy for 51 and 6lm [123]. Of crucial importance is the observation that the 
position of the third peak appears to remain insensitive to other cosmological quantities. 
This is because at small scales, well inside the horizon, ^ is usually negligible, hence we 
recover the unperturbed harmonic oscillator equation for Eq. H4.5|) that describes oscil- 
lations with zero point offset. Hence we can make use of this fact to test dark energy 
models [124,125]. As we shall discuss in more detail in Chapter |S1 the quintessence field 
can leave a distinctive signature on the shape of the CMB power spectrum through both 
the early integrated Sachs- Wolfe effect (ISW) and the late one [126]. The former is impor- 
tant if the dark energy contribution at the last scattering surface is not negligible [100, 118] 
or in non-minimally coupled models [82, 127, 128], whereas the late ISW is the only effect 
in models with Q,q ~ such as those described by our parameterized quintessence poten- 
tial [129]. However the late ISW produces an imprint on the CMB power spectrum of the 
order of 10% and therefore is not detectable with the pre-WMAP measurements. In such 
a case an accurate determination of the position of the Doppler peaks is more sensitive to 
the actual amount of dark energy. 

4.4 Likelihood analysis and results 
4.4.1 Constraints from supernovae 

We want to constrain the set of parameters a, (3 and VIq confined in the range: a G (1, 10), 
(3 G (0, 10) and VIq G (0, 1), subject to the assumption of a flat universe. This choice of 
the parameters allows us to account for a large number of models. We use the Sn la data 
fit C of Perlmutter et al. (1999) [53], that excludes 4 high redshift data points, which are 
very likely reddened by their host galaxies. The magnitude-redshift relation is given by: 



where ^ is the 'Hubble constant free' absolute magnitude and Dl{z) = Hodiiz) is the 
free-Hubble constant luminosity distance. In a flat universe 



where tq is the conformal time today and t{z) is the conformal time at the red-shift z 
of the observed supernova. Both of these quantities are calculated solving numerically 



m{z) = 5 log Dl{z, a, (3,^q) + 



(4.8) 



dUz) = {to - t{z)){1 + z) 



(4.9) 
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Figure 4.4: Fractional Quintessence energy density likelilioods, (a) for Sn la, (b) for the 
combined CMB peaks and (c) for the combined data sets. 

Eq. (|2.8() and Eq. (|2.9() for each value of a, (3 and Og. In ^ we neglect the dependence 
on a fifth parameter (a in [53]) and assume it to be 0.6, the best fit value in [53]. We then 
obtain a gaussian likelihood function ^"^"(a, /3, Jig), by marginalizing over In figure 
14.4b we present the one-dimensional likelihood function normalized to its maximum value 
for r^Q. When considering only Sn la data, there is a maximum at fig = 1, such a high 
value of Qq is required if we impose the constraint 0^ = and it is in agreement with the 
analysis in [130]. In figure l^3k we present the likelihood contours in the a — (3 parameter 
space, obtained after marginalizing over Qq. Note that all values are allowed at the 2a 
level. The confidence regions for the Sn la data correspond to Quintessence models with 
Wq < —0.4 for Qq = 0.6, an upper limit that agrees with those found in [53, 114]. 

4.4.2 Constraints from Doppler peaks and Sn la 

We now compute the position of the three Doppler peaks li, h and Z3 using Eq. 1)4.7(1 . In 
addition to the parameter space used in the supernovae analysis we consider the physical 
baryon density and the scalar spectral index varying respectively in the range Q.},h? G 
(0.018,0.026) and n G (0.9,1.1). The Hubble constant is set to h = 0.70 in agreement 
with the recent HST observations [29] . It is worth remarking that the baryons density and 
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Figure 4.5: Likelihood contour plots for Sn la, I, II and III acoustic peaks. The blue 
region is the 68% confidence region while the 90% is the light blue one. For the Sn la 
the white region correspond to 2a. The position of the third CMB acoustic peak strongly 
constrains the acceptable parameter space. 

Hq are degenerate with the quintessence parameters. In fact increasing ri^/i^ reduces the 
value of Cs causing a shift of the CMB peaks towards large multipoles. The same effect 
occurs for low values of the Hubble constant, therefore the results of this analysis will 
depend on the HST prior. The predicted peak multipoles are then compared with those 
measured in the Boomerang and DASI spectra [131]. Note, that the third peak has been 
detected in the Boomerang data but not in the DASI data. In table 14.11 we report the 
position of the peaks with Icr errors estimated from the Boomerang, DASI and WMAP 
data with a model independent analysis. Since the errors are non Gaussian (see also [132]), 
to be conservative we take our Icr errors on the data to be larger than those reported 
in [131], so that our analysis is significant up to 2a. We then evaluate a gaussian likelihood 
function ^^^"'^'^{a, /3, Oq, Vli,h? ,n). The combined one-dimensional likelihood function for 
the peaks is shown in figure IT^ . where we find Jig = CGOitQ^Q. The likelihood for all 
the data sets combined is shown in figure where we find Q.q = 0.75ibQQg. These 

results are in agreement with the analysis in [21, 133, 148]. The likelihood functions. 
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Peak Boomerang DASI 



WMAP 



h 



213 



+10 

'-13 



202 ± 15 220.1 ±0.8 




548 ± 10 



546 ± 10 




Table 4.1: Location of the CMB peaks inferred from the power spectrum of Boomerang 
and MAXIMA data using a model independent procedure (from [131]). 

combining all the CMB peaks data, for the scalar spectral index and the physical baryon 
density are shown in figure ITBl Since the dependence of the peak multipoles on rj^/i^ and 
n is weak, it is not possible to obtain some significant constraints on these cosmological 
parameters. In figure l4.5bM.5l i we plot the two-dimensional likelihood function in the 
plane a — (3 for each peak, obtained after having marginalized over Og, O^/i^ and n. Their 
shape reflects the accuracy in the estimation of the position of the peaks. Actually the 
first one is very well resolved, while we are less confident with the location of the second 
and third peak. Therefore their likelihoods are more spread and flat in the a — (3 plane. 
The 1(7 confidence contour (figure H3b ) for the first acoustic peak constrains the slopes 
of our potential in the range: 3 < q < 10 and 1 < /3 < 3. In particular the likelihood 
has a maximum at a = 9 and (3 = 2, corresponding to an equation of state Wq = — 0.8 
for Qq = 0.7, in agreement with the recent analysis in [133]. However, the second and 
third peaks constrain a complementary region where the equation of state is compatible 
with the cosmological constant value. Therefore the effect of including all the data in the 
likelihood analysis is to move the constraint from models with Wq ~ —0.8 to models with 
an equation of state Wq ~ — 1. The constraints on the slopes of the quintessence potential 
allows us to infer an upper limit on the present value of the equation of state only after 
having fixed the value of Qq. In fact the scalar field potential Eq. 1)4.4(1 is fully specified 
by three numbers that are M, a, [5. However the map (M, a, j3) — > nfq is not one-to-one 
and therefore it is not possible to transform the likelihood ^(a,/3, ilg) ^{iJq,Q.q). 
Therefore there is a substantial difference between inferring the bounds on the dark energy 
equation of state by directly constraining its present value and constraining the scalar field 
potential. As we can see in figure 14.81 the values of a and /?, allowed by the likelihood 
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Figure 4.6: One-dimensional likelihood for n and Qfo/i^. 

including all the data (figure HTj) . correspond to our models with values of Wq in the range 
-1 < < 0.93 at 2(7 for our prior probability ^Iq = 0.75. Assuming smaller values of 
Qq relaxes the constraints on Wq, for instance for Qq = 0.5 we have —1 < Wq < —0.6 at 
2(7. The results of our analysis can be interpreted as follows. At the la level the position 
of the first peak is inconsistent with the position of the other two. A possible explanation 
of this discrepancy is that the multipoles I2 and I3 are less sensitive to small shifts induced 
by the dependence on Ofe/i^ and 7i. Therefore we can obtain a different constraint on the 
dark energy equation of state if we consider the peaks individually. It is worth remarking 
that the position of the first peak in the Boomerang data prefers slightly spatially closed 
models. Having assumed a flat geometry affects our upper bounds on the slope of the 
potential in a region corresponding to Wq > — 1. On the other hand the location of the 
third peak inferred in [131], is at l-^ = 845ib23, for which values of Wq ~ — 1 fit this 
multipole better than models with Wq > — 1. In fact the peaks are shifted toward larger 
multipoles as Wq approaches the cosmological constant value. This is because, in models 
with Wq ~ — 1 the Universe starts to accelerate earlier than in those with Wq > — 1, 
consequently the distance to the last scattering surface is farther and hence the sound 
horizon at the decoupling is projected onto smaller angular scales. A possible source of 
bias in our analysis is due to the prior on the value of f^o, in fact assuming h = 0.70 
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Figure 4.7: Two-dimensional likelihood for Sn la and CMB with 1 (dark blue) and 2a 
(light blue) contours. 

implies that models with a large negative equation of state and large value of are 
preferred. However our results are consistent with other studies which assume similar 
priors. For instance in [133], using the complete set of available CMB measurements the 
slope of the Inverse Power Law potential has been constrained to be a = CSiLqI under 
the prior h = 0.65, corresponding to Wq = —0.82±Q'\f. Similarly in [134], by making use 
of the CMB, Sn la and large scale structure data, it has been found wq < —0.85 for a 
constant equation of state. A time dependent dark energy equation of state, characterized 
by a late time transition, has been constrained using all cosmological data in [135], where 
it has been found Wq = — l^'^'^. 

4.5 Discussion 

The analysis of Sn la data and the location of the CMB peaks in Boomerang and DASI 
data constrain the slope of a general quintessence potential in a range of values such that 
the quintessence field is undergoing small damped oscillations around a minimum of the 
potential (see figure Such a situation implies that the value of the quintessence 

equation of state has to be very close to —1. However this does not imply that the dark 
energy is the cosmological constant. What our result means is that the quintessence 
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Figure 4.8: Contour values of in the a — [5 plane for = 0.5 (upper panel) and 
VLq = 0.75 (upper panel). 

field today is evolving in a nearly vacuum state. In fact models where the scalar field is 
rolling down a very flat region of the potential can fit the data and pass the constraints 
even though they are not included in our analysis. Another important caveat is that this 
study does not take into account quintessence scenarios where the contribution of the dark 
energy is not negligible. In such a case we would have to take into account physical effects 
not only on distance measurements, but on the structure formation process itself (see for 
instance [104,136]). These effects however are not detectable with the level of accuracy 
provided by the recent balloon and ground CMB experiments. A revised analysis of the 
Boomerang data has been recently published in [22], the position of the peaks has been 
determined with a better accuracy, in particular the third peak is shifted at a smaller 
multipole than previously detected. This slightly changes our results, relaxing the lower 
bounds on the slope /3. Here we want to stress the fact that it seems likely the present 
value of the dark energy equation of state is close to —1. However the possibility that the 
equation of state was largely different from —1 cannot be excluded, as it is shown by the 
class of models that best fit the Sn la data and the CMB peaks. The new generation of 
satellite experiments we will probably be able to detect the specific effects produced by a 
time varying dark energy component. 



Chapter 5 

A model independent approach to 
the dark energy equation of state 



The detection of time variation in the dark energy equation of state can be considered 
as the 'smoking gun' for the cosmological constant scenario. The new generation of ex- 
periments in cosmology wiU provide high precision measurements that in principle can 
discriminate between different dark energy candidates. On the other hand it is unrealistic 
to assume we can infer some general information constraining some particular quintessence 
model. In fact a plethora of scalar field potentials have been proposed, all leading to sim- 
ilar late time behaviour of the universe, hence to a lack of predictability. In this Chapter, 
we first review some of the methods that have been proposed in the literature to constrain 
the time features of a general quintessence component. We then introduce a model inde- 
pendent approach which simply involves parameterizing the dark energy equation of state 
in terms of known observables [137]. This allows us to analyse the impact dark energy has 
had on cosmology without the need to refer to particular scalar field models and opens up 
the possibility that future experiments will be able to constrain the dark energy equation 
of state in a model independent manner. 



5.1 The effective equation of state 

Measuring the present value of the dark energy equation of state can distinguish between 
a cosmological constant scenario and a quintessential component if fWg > — 1. Using the 
tracker properties of the quintessence field, the most reliable method to infer a bound 
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(A4+"/Q")et«' [96] A4(coshAQ- 1)P [138] 

^4(ga^Q ^ g/3KQ) [;^oO] K'^e-^^{l+As{nvQ) [139] 

A4[(Q-5)" + A]e-^« [118] A4[1 + cos(Q//)] [74] 

^'(^ + ^) A4(eVQ-l) [79] 

Table 5.1: Quintessence potentials that have been used in the litterature. 

on ioq would be to constrain directly the quintessence potential. However this could 
result in a difficult challenge, for instance in Table we show a list of proposed tracker 
quintessence potentials. Therefore, because there are no fundamental physical principles 
that can specify the nature of the dark energy, we are left with a potentially infinite set of 
families of theoretical models to compare with the data. As a consequence of this it is to 
be hoped that a more general way of constraining the dark energy can be developed. A 
simple method to take into account a large number of dark energy models is to consider 
a constant effective equation of state defined by Eq. 1)4. 3|) . Using the energy conservation 
equation, the redshift evolution of the dark energy density is given by 

pUz) = p^,,{l + zf^^+^^ff\ (5.1) 

The effective equation of state appears through the Hubble equation in the coordinate 
distance defined by: 

1 dz' 

Hence the value of Weff can be constrained by measurements of the luminosity distance 
and the angular diameter distance, that in a fiat space-time read as di{z) = (1 + z)r{z) 
and d^iz) = r{z)/{l + z) respectively. There is a general consensus that an accurate 
determination of from high redshift Sn la measurements with the proposed SNAP 
satellite [140] can determine Weff to better than 5 per cent [141]. In such a case if "u^e// 
is much larger than —1 it will be possible to rule out the ACDM model. However if 
Weff ~ — 1 the results of the data analysis have to be carefully considered. In fact it has 
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been shown by Maor et al. [142] that assuming a constant t^e// introduces a bias that 
leads to misleading conclusions about the properties of the dark energy. In order to make 
this point clear, let us consider a sample of mock data generated from a dark energy model 
with an equation of state of the form 

wq{z) = wq + wiz, (5.3) 

with wq < —1/3 and wi > 0. The constraints on a constant equation of state w^ff 
obtained from the statistical analysis of such a sample, will push the bounds towards values 
of Weff < Wq. This effect has been noticed also in other works, see for instance [145,146] 
and can be explained as follows. For a model with wq given by Eq. ()5.3() the energy 
density evolves as: 

As we may note, comparing Eq. ()5.1() with Eq. 1)5. 4|) . if t^o < and wi > 0, fitting the data 
with a constant equation of state roughly requires Weff ^ wq — wi < wq. In other words 
the fitting procedure will pick out values of wq that are more negative than the fiducial 
one. Hence a time varying dark energy component can be confused with a cosmological 
constant, while allowing Weff to vary in a range of values more negative than -1, will 
shift the best fit models towards Weff < — 1. Therefore it is very possible that the claims 
for the existence of phantom matter (with wq < —1) [143,144] are just a consequence of 
assuming a constant equation of state in the data analysis. 



3wiz 
l + z 



(5.4) 



5.2 Cosmological distance fitting functions 

Several methods to constrain the time dependence of the dark energy equation of state have 
been discussed in the literature [115, 145-153]. Since the dark energy became dominant 
only recently, it is reasonable to concentrate on local redshift measurements, such as 
the Sn la. The coordinate distance and the Hubble parameter are uniquely related by 
H{z) = {dr / dz)^^ . Thus we can estimate r{z) from cosmological distance indicators and 
unambiguously calculate H{z) and dH/dz. This allows us to reconstruct w{z) provided 
the value of Vim is given. In fact the Hubble equations in the presence of matter and a 
scalar field are: 



H 



2 



(5.5) 
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and 

H = -A7TG{pr,^+Q^), (5.6) 

where the dot is the derivative with respect to the time. These equations can be rewritten 
in the form: 

mf^''^ ~ Hi ml \dx) 2^^™"^ ' ^^-^^ 



SvrG fdQy 2 O^x 



?,Hl \dx J mix ' 
where x = 1 + z. Using the definition of the scalar field equation of state we find: 



2x dlnH -j^ 

w{z) = ^-^1 . (5.9) 

^0 

The idea is then to introduce a fitting function for the coordinate distance, so that once its 
parameters have been determined from the data analysis, its first and second derivatives 
with respect to the redshift can be analytically calculated. The authors of [115] suggest 



-Hi 



ay/x — 1 + a 



(5.10) 



Px + 'jy/x + 2 — a — /? — 7 
where a, [5 and 7 are fitting parameters, while Huterer & Turner [147] consider both a 
polynomial fitting function 

rf'\z) = ^CiZ^ (5.11) 
i 

and a Pade' approximate 

r/-(z) = -4^- (5.12) 



It has been pointed out in [146] that such a general method fails to accurately reproduce 
the correct behaviour of given scalar field models, even assuming strong priors on the value 
of $7^. The difficulty arises because the formula Eq. (|5.9|) depends on the derivatives of 
r[z). Therefore even though r^'^^{z) can mimic the r{z) of a specific model, this may be 
not true for its derivatives. 
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5.3 Statefinder method 

A different approach has been developed by Sahni et al. [154,155]. They suggest a method 
to distinguish amongst different dark energy scenarios by measuring geometrical parame- 
ters that are constructed in such a way they are more sensitive to the expansion rate of 
the Universe. Such parameters, called Statefinder pair, are defined by: 

a r — 1 , , 

,s = - — — . (5.13) 



am' 3(g-l/2)' 

For a flat Universe the relation to the dark energy equation of state w and its time 
derivative w is: 

9 ?> w 

r = l + -Vldew{l + w)--nde—, (5.14) 

1 w , ^ 

s = l + w-- — , 5.15 
6 wH 

where H is the Hubble parameter. In general the value of such parameters evolves with the 
time, but for the cosmological constant {w = —1) we have (r, s) = (1,0) at all the times. 
Moreover, since (r, s) are geometrical quantities, they can naturally take into account 
braneworld models in which the accelerated expansion is caused by five dimensional gravity 
effects and which lead to modifications of the Hubble equations. The Statefinder pair can 
be estimated using the following ansatz for the Hubble parameter, 

H{x) = Ho [nmx^ + Ai+ A2X + A^x"^] ^ , (5.16) 

where x = 1 + z, and Ai + A2 + A^ = 1 — r^^- In such a case (r, s) are of the form: 

^ ^ Al + n^x^ ,^ 

^ VLmX''^ + Ai+ A2X + A^.x'^' 

2{A2X + Asx^) 

3{Mi+2A2x + A3X^)' ^ ' 

The value of the fitting parameters Ai can be inferred from a likelihood analysis of cosmo- 
logical distance measurements, such as the Sn la. From Eq. (|5.17|1 it appears obvious that 
an accurate estimation of r needs a precise knowledge of the amount of non-relativistic 
matter. In [155] it has been shown that imposing a Gaussian prior probability distribution 
on Qrn. with variance (Tq^^ = 0.05, the high redshift Sn la measurements provided by the 
proposed SNAP satellite will accurately determine the Statefinder pair offering the chance 
to distinguish between different dark energy models. Therefore the Statefinder method 
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seems to be highly sensitive to the nature of dark energy and avoids the main difficul- 
ties of the reconstructing procedure. In fact it only needs to constrain fitting parameters 
that appear in the Hubble equation and does not constrain directly the equation of state 
that affects the expansion rate apart from through a time integrated effect. A potential 
limitation of such an approach is that the determined values of (r, s) do not give us any 
information about the physical properties of the dark energy, unless their values are a 
priori known for specific classes of models. 

5.4 Low redshift parameterization 

In the absence of theoretically motivated dark energy models we can consider a low red- 
shift expansion for 'Wq{z) and the parameters of such an expansion can be interpreted as 
fundamental dark energy parameters. Therefore the constraints on a parameterization of 
the equation of state would provide a model independent determination of the dark energy 
properties. Such a parameterization should satisfy the following requirements: 

• it should depend on a minimal number of parameters 9i that specify the physical 
properties of the dark energy; 

• the functional form of the parameterization has to be such that for a given set {6i) 
the behavior of the equation of state wq{z,9) reproduces the w{z) predicted by 
proposed quintessence models and can account also for more general cases. 

In [146, 148] it was suggested the use of a logarithmic expansion, 

Wq{z) = Wo - a log{l + z), (5.19) 
while the authors of [145, 151, 153] used a polynomial fitting function, 

w'PQ = ^Ciz\ (5.20) 

i 

The general drawback of using a redshift expansion formula for the dark energy equation 
of state is that the number of fitting parameters is now larger than assuming a simple 
time constant behaviours and consequently the degeneracy in the dark energy parameter 
space is enlarged. In fact the coordinate distance is related by a multi-integral expression 
to the equation of state and therefore widely different w{z,6i) can have the same multi- 
integral value. There is a general consensus that assuming strong priors on the value of 
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Figure 5.1: Low redshift approximation to the quitenssence equation of state for different 
types of potentials: Inverse Power Law (red), Inverse Power Law times Exponential (blue) 
and Two Exponential (magenta). The solid line corresponds to the exact numerical so- 
lution of the Klein-Gordon equation, while the dash line is the approximated behaviour 
provided by Eq.^J^. 

i^rn it will be possible to tightly constrain wq = 8i, but it will be more difficult to put 
tight bounds on the other parameters. The solution would be to combine different data 
sets that can break the degeneracy in the parameter space. The expansions specified by 
Eq. (|5.19|) and Eq. (|5.2U() both have limited applicability. For instance Eq. ()5.19() can take 
into account only for quintessence models with slowly varying behaviours and breaks down 
at large redshifts, as we can see in figure l5!Tl On the other hand the polynomial expansion 
Eq. (|5.20() introduces a number of unphysical parameters whose values are not directly 
related to the properties of a dark energy fluid. The consequence is that, unless we pay 
the cost of considering very unphysical values, their application is limited to low redshift 
measurements and cannot be extended to the analysis of high redshift data sets such as 
the CMB data. An interesting alternative to the fitting expansion approach, has recently 
been proposed in [156]. The authors of this study developed a method to reconstruct the 
time behavior of the equation of state from cosmological distance measurements without 
assuming the form of its parameterization. In spite of the efficiency of such an approach. 
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it does not take into account the effects of the possible clustering properties of dark energy 
which become manifest at higher redshifts. Hence its application has to be limited to the 
effects dark energy can produce on the expansion rate of the universe at low redshifts. On 
the other hand, it has been argued that dark energy does not leave a detectable imprint 
at higher redshifts, since it has only recently become the dominant component of the 
universe. Such a statement, however, is model dependent, on the face of it there is no 
reason why the dark energy should be negligible deep in the matter dominated era. For 
instance CMB observations constrain the dark energy density at decoupling to be less then 
10 per cent of the critical one [104]. Such a non negligible contribution can be realized in 
a large class of models and therefore cannot be a priori excluded. Consequently it is of 
crucial importance to find an appropriate parameterization for the dark energy equation 
of state that allows us to take into account the full impact dark energy has on different 
types of cosmological observations. 

5.5 An exact parameterization for the dark energy equation 
of state 

Our goal is to determine an appropriate analytical form of the equation of state WQ{a) 
valid at all redshifts in terms of physical quantities, so that it can describe a general 
fluid and take into account most of the proposed dark energy models. In Chapter 2 we 
have seen that the specific evolution of 'WQ{a), depends on the shape of the potential, 
however there are some common features in its behaviour that can be described in a 
model independent manner and which allow us to introduce some physical parameters. 
As we have seen in Chapter 2 a large number of quintessence models are characterized by 
the existence of the 'tracker regime'. It consists of a period during which the scalar field, 
while it is approaching a late time attractor, evolves with an almost constant equation of 
state whose value can track that of the background component. Here we consider a broad 
class of tracking potentials. These include models for which WQ{a) evolves smoothly, 
as with the inverse power law [73], V{Q) ^ (INV) and the supergravity inspired 

potential [96], V{Q) ~ l/Q^e^'/^ (SUGRA). Late time rapidly varying equation of states 
arise in potentials with two exponential functions [100], V ~ e~"*^ + e^^ (2EXP), in the 
so called 'Albrecht &: Skordis' model [118] (AS) and in the model proposed by Copeland 
et al. [101] (CNR). To show this in more detail, in figure EH we plot the equation of state 
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Figure 5.2: Evolution of wq against tiie scale factor for an Inverse Power Law model (solid 
blue line), SUGRA model ( [96]) (dash red line), two exponential potential model ( [100]) 
(solid magenta line), AS model ( [118]) (solid green line) and CNR model ( [101]) (dot 
orange line). 

obtained by solving numerically for each of these potentials the scalar field equation of 
motion defined by Eq. and Eq. 

There are some generic features that appear to be present, and which we can make 
use of in our attempts to parameterize wq. For a large range of initial conditions of the 
quintessence field, the tracking phase starts before matter-radiation equality. In such a 
scenario WQ{a) has three distinct phases, separated by two 'phase transitions'. Deep into 
both the radiation and matter dominated eras the equation of state, WQ{a), takes on the 
values Wq and Wq respectively, values that are related to the equation of state of the 
background component wb through: 



where F = V"V/{V')'^ and V = dV/dQ etc. For the case of an exponential potential, 
F = 1, with Wq = Wb, but in general wq 7^ wb- Therefore if we do not specify the 
quintessence potential the values of Wq and Wq should be considered as free parameters. 
The two transition phases can each be described by two parameters; the value of the scale 
factor ac™" when the equation of state wq begins to change and the width A,.,^ of the 
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Figure 5.3: Plot of WqIo) fEq. 15.221) best fit for different potentials. 

transition. Since T is constant or slowly varying during the tracker regime, the transition 
from Wq to Wq is always smooth and is similar for all the models (see figure 1^^ . To be 
more precise, we have found that aJJ ~ 10~^ and ~ 10~^ during this transition, the 
former number expected from the time of matter-radiation equality and the latter from 
the transition period from radiation to matter domination. However, when considering 
the transition in wq from to the present day value nfq , we see from figure 15.21 that 
this can be slow (0 < / /S.m < 1) or rapid (a™/Am > 1) according to the slope of the 
quintessence potential. For instance in models with a steep slope followed by a flat region or 
by a minimum, as in the case of the two exponentials, the AS potential or the CNR model, 
the scalar field evolves towards a solution that approaches the late time attractor, finally 
deviating from the tracking regime with the parameter F rapidly varying. In contrast the 
inverse power law potential always has a slower transition since F is constant for all times. 
Given these general features we conclude that the behavior of WQ{a) can be explained in 
terms of functions, WQ{a), involving the following parameters: Wq, Wq, Wq, a™ and A^. 
The authors of [135] have recently used an expansion in terms of a Fermi-Dirac function 
in order to constrain a class of dark energy models with rapid late time transitions. 

In what follows we find that a generalisation of this involving a linear combination of 
such functions allows for a wider range of models to be investigated. To be more precise. 
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.0 w^Q a™ A 



m 



INV -0.40 -0.27 -0.02 0.18 0.5 

SUGRA -0.82 -0.18 0.10 0.1 0.7 

2EXP -1 0.01 0.31 0.19 0.043 

AS -0.96 -0.01 0.31 0.53 0.13 

CNR -1.0 0.1 0.32 0.15 0.016 

Table 5.2: Best fit values of the parameters of the expansion ()5.22|1 . 

we propose the following formula for ^^(a): 

w^Qia) = FiU{a) + F2fm{a) + F3, (5.22) 

with 

fr,m{0') — ^_^r,7?i • (5.23) 

The coefficients Fi, F2 and F3 are determined by demanding that WQ{a) takes on the 
respective values Wq, Wq, Wq during radiation (ur) and matter (om) domination as well 
as today (ao)- Solving the algebraic equations that follow we have: 

i^Q - '^W (/m(ao) - fmiar)) " jw'^Q " W^q) (/m(am) " fmjar)) 
ifriam) - friar)) ifm{ao) " fm{ar)) " (/r(ao) " friar)) {fm{am) " fm{ar)) ' 

(5.24) 



fm{ao) — fm{ar) /m(oo) — fmiar) ' 

F3 = w'q- Fi friar) - F2 fmiar), (5.26) 

where qq = 1, and the value of and a,, and be arbitrarily chosen in the radiation and 

matter era because of the almost constant nature of wq during those eras. For example 
in our simulations we assumed = 10~^ and = 10~^. In table we present the 
best fit parameters obtained by minimizing a chi-square for the models we have considered 
and in figure we plot the associated fitting functions WQia). It is encouraging to see 
how accurately the Fermi-Dirac functions mimic the exact time behavior of wqia) for the 
majority of the potentials. 

In figure Em we plot the absolute value of the difference Aw(a) between wqia) and 
Wqia). The discrepancy is less then 1% for redshifts z < 10 where the energy density of the 
dark energy can produce observable effects in these class of models and it remains below 9% 
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Figure 5.4: Absolute value of the difference between WQ{a) and WQ{a) for the models of 
figure 15.21 

between decoupling and matter-radiation equality. Only the CNR case is not accurately 
described by ^^(a) due to the high frequency oscillations of the scalar field which occur at 
low redshift as it fluctuates around the minimum of its potential. In fact these oscillations 
are not detectable, rather it is the time-average of WQ{a) which is seen in the cosmological 
observables, and can be described by the corresponding ^^(a). There are a number of 
impressive features that can be associated with the use of WQ{a) in Eq. 1)5.22(1 . For instance 
it can mimic the behavior of more general models. As an example of this in figure 1^31 we 
plot WQ{a) that approximately corresponds to three cases: a K-essence model [157] (blue 
solid line); a rapid late time transition [158] (red dash-dot line) and finally one with an 
equation of state Wq < — 1 (green dash line). The observational constraints on Wq, Wq , 
Wq, a™ and Am lead to constraints on a large number of dark energy models, but at the 
same time it provides us with model independent information on the evolution of the dark 
energy. It could be argued that the five dimensional parameter space we have introduced 
is too large to be reliably constrained. Fortunately this can be further reduced without 
losing any of the essential details arising out of tracker solutions in these Quintessence 
models. In fact nucleosynthesis places tight constraints on the allowed energy density 
of any dark energy component, generally forcing them to be negligible in the radiation 
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era [103,159]. The real impact of dark energy occurs after matter-radiation equality, so 
we can set Wq = Wq. Consequently we end up with four parameters: Wq, Wq, and 
Am- Although they increase the aheady large parameter space of cosmology, they are 
necessary if we are to answer fundamental questions about the nature of the dark energy. 
The parameters make sense, if WQ{a) evolves in time, we need to know when it changed 
(a^), how rapidly (Am) and what its value was when it changed (wq). Neglecting the 
effects during the radiation dominated era it proves useful to provide a shorter version of 
Eq. 1)5.22(1 . in fact since we can neglect the transition from radiation to matter dominated 
eras, then the linear combination Eq. (|5.22j) can be rewritten as ^: 

771 

1 _|_ p 1 P "^771 

^^(«) = ^'q + K - ^'q) X \^ X (5-27) 

1 _l_ e A,„ 1 — e ^™ 

As we can see in figure EEl the relative difference between the exact solution WQ{a) of the 

Klein-Gordon equation and Eq. I|5.27() is smaller than 5% for redshifts z < 1000, therefore 

it provides a very good approximation for the evolution of the quintessence equation of 

state. Both Eq. ()5.22() and Eq. ()5.27() are very useful in that they allow us to take into 

account the clustering properties of dark energy (see for instance [160]) and to combine 

low redshift measurements with large scale structure and CMB data. 



^We thank Eric Linder for pointing this out to us. 
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Figure 5.5: Time evolution of WQ{a) as in the case of K-essence (blue solid line), late time 
transition (red dash-dot line) end with Wq < — 1 (green dash line). 




Figure 5.6: Absolute value of the difference between tcg(a) and the low redshift formula 
Eq. ()5.27|) for the models of figure [2121 



Chapter 6 

Dark energy effects in the Cosmic 
Microwave Background Radiation 



The physics of the Cosmic Microwave Background radiation has been deeply studied during 
the past 30 years. The theoretical developments in understanding the different processes 
responsible for the temperature anisotropies have given to us the possibility to accurately 
predict the spectrum of the CMB fluctuations for a given cosmological model. The exis- 
tence of these cosmological temperature fluctuations was initially confirmed by the COBE 
satellite [163], but it is only with the recent generation of balloon and ground experiments, 
such as Boomerang [19], Maxima [23] and DASI [24], that the observations reached the 
level of accuracy necessary for testing different cosmological scenarios. As we have reviewed 
in Chapter ^ these measurements have detected the multiple peak structure of the CMB 
power spectrum, providing an accurate determination of the cosmological parameters. The 
improvements in the CMB experiments are mainly due to the high performance of the new 
high frequency microwave detectors. The recent WMAP satellite, using this technology 
has started a new generation of satellite experiments that will measure the temperature 
and polarization anisotropies close to the theoretical 'cosmic variance' limit. It is there- 
fore of crucial importance to study the imprint dark energy leaves in the anisotropy power 
spectrum. In this Chapter we start reviewing the main concepts of the CMB physics. 
Using a model independent approach we study the impact of different dark energy models 
in the CMB spectrum. In particular we find that features of the dark energy equation 
of state can leave a characteristic imprint in the Integrated Sachs- Wolfe effect [161,162]. 
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We will identify the dark energy models that are distinguishable from the cosmological 
constant. In conclusion we will show the limits on a class of dark energy models that can 
be obtained by cosmic variance limited experiments. 



6.1 A beginner's guide to CMB physics 

The starting point of the CMB anisotropy calculation is the kinetic theory of photons in 
a perturbed space-time. Here we shall briefly review the main concept of the subject and 
for more authoritative reviews we refer to [121,126,164]. 

6.1.1 Basic equations 

Let us consider a flat pertubed background. In the Newtonian gauge the metric element 
is 

= + 2^)dr]^ + (1 - 2<^)dx'dx^], (6.1) 

where ^ and $ are the metric perturbations and rj is the conformal time. The properties of 
a gas of photons propagating in this space-time are described by the distribution function 
/(x*, Pi) that depends on the three spatial coordinates x* and the three conjugate momenta 
Pi. In the Newtonian gauge the relation between Pi and the proper momentum pi measured 
by an observer at a fixed spatial coordinate is given by Pi = a(l — ^)pi. In order to 
eliminate the metric perturbations from the definition of the momenta it is convenient 
to replace Pi with qi = api, where qi = q^i, q and -ji being the module and the cosine 
directions respectively. In the early Universe the photons propagate in a hot ionized 
medium and interact through the Compton scattering with free electrons and ions. As 
a consequence of this the distribution function is conserved along the perturbed geodesic 
apart for a coUisional term. The evolution of the distribution function is described by the 
Boltzmann equation: 

^ ^ ^ + + + = c[f] (6 2) 

dr] df] dx^ drj dq drj dj^ dr] ' 

where C[f] accounts for the Compton scattering and the geodesic equation 

Idq _^ ^d'^ 
q dr] * dxi 

In addition the linearized Einstein equations determine the equations for the metric per- 
turbations and ^ in term of the perturbations in the multiple fluids system. At the 
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zeroth order we can consider the Universe as being unperturbed, in sucli a case / is inde- 
pendent of and 7^ and if tlie collision between photons and charged particles conserve 
energy, then q = ap is time independent and consequently / is as well. This is a good 
approximation, since the free electrons are non-relativistic and the Compton scattering 
is primarily Thomson scattering in which momentum and not energy is transferred. The 
equilibrium distribution function / is 



where the factor 2 accounts for the number of polarization states of the photon and Tq is 
the temperature of the radiation today. The effects of the perturbed background can be 
described by expanding the distribution function / about the equilibrium state, 

f = fo{q) + fi{v,^\li,q), (6.5) 

where fo{q) is given by Eq. (|6.4() and /i is the perturbed part of the distribution function. 
Substituting Eq. (|6.5j) into Eq. H().2j) and taking only the first order terms we obtain: 

drj ~^ ^ dx^ ~^ dq ^ ^ dx^ ^ [/] ( ) 

For practical purposes it is convenient to write the perturbed part of the distribution 
function in term of the brightness function 0(?/, x*,7j), that is the fractional perturbation 
in the effective temperature of the photons, T = ro(l + 0). Expanding / around To(l + 0) 
we have: 

/ 



^^P [T(,(i+e)] 



2 „ q exp(^) 



exp (^) - 1 



2 2 „^°^ 



exp(^)-l 



Tf 







= /o - gf^e. (6.7) 

aq 

Comparing Eq. (|6.5|) with Eq. 1)6. 7|) we obtain 

fl = -qdfo/dqe. (6.8) 

Substituting Eq. 1)6. 8|) into Eq. 1)6. 6() and taking the Fourier transform the Boltzmann 
equation becomes [18]: 

G + iki2{e + '^) = -^ + Sc, (6.9) 
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where /i = and 



(6.10) 



where f = XeUeaTCL is the differential optical depth, with the ionization fraction, ng the 
electron number density and or the Thomson cross section; ©o is the isotropic component 
of G, V^* is the baryon velocity and Hr is the anisotropic stress perturbation for the 
photons. Hence the collisional term couples the evolution of the photon perturbation to 
that of the baryons. The equations for the baryon perturbations are obtained by linearizing 
the conservation equation of the stress energy tensor: 



4 = -HVb - Vr) - 3^, 

Vb = --Vb + + f{Vr - Vb)/R, 



(6.11) 
(6.12) 



where Vr is the velocity perturbation of the photons. It is useful to expand the brightness 
function in a Legendre series: 



& = J2^-iy{2i+i)@iPiif^). 



(6.13) 



1=0 



Using the definition of the photon stress energy tensor in term of the distribution function 
/, it can be shown that ©o = 5^/4 and ©i = Vr- Substituting Eq. ()6.13|) in Eq. 1)6. 9 j) and 
using the explicit form of the Compton scattering term Eq. ()6.1U|1 we obtain the hierarchy 
equations [165,166]: 



©n 



©1 = A; 



eo = k 



e, = k 



©0 + ^ - -©2 
5 



t(0i - Vb), 



9 



21 - 1 



-© 



i-i 



10 

1±1 
21 + 3 



f©2, 



© 



/+1 



fQi (/ > 2). 



The anisotropy power spectrum is defined as: 



2/ + 1 ^ J_ 
47r ' ~ 27r2 



drj 



dkk^\@l{k,r])\'^ 



(6.14) 
(6.15) 

(6.16) 

(6.17) 

(6.18) 



k 21 + 1 

From Eq. (|6.9j) we note there are two sources of anisotropy formation, the Compton 
scattering that couples the photons to the baryons and the gravitational effect produced 
by the presence of density fluctuations in all the matter components. 
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6.1.2 CMB anisotropies 
Acoustic Oscillations 

Before recombination the differential optical depth f is very large and the scattering 
between photons and baryons is extremely rapid and efficient. In this regime, called 
the tight coupling limit {f/k >> 1), the photons and the baryons behave as a single 
fluid. Because of this the baryons and photons velocity are the same at the zeroth order, 
Vh ^ Vr = &!■ In such a case Eq. H6.12() becomes: 

ei = -|ei + fc* + ^(ei-i4), (6.19) 

where R = 3ph/4:pr is the baryon to photon ratio. Inverting this equation in terms of 
f(0i — Vf,) and substituting in Eq. (|6.15|) we obtain the iterative first order solution, 

60 = -|e, - *, (6.20) 

These equations can be combined into a single second order equation [167], 

®o + -T^Qo + = ^(^)' (6.22) 

al + R 

where = ^{l + R)~^. The source function F{r]) arises from the gravitational potentials 
and is given by 

al + R 6 

From Eq. (|6.22() we can distinguish three different contributions to the evolution of the 
isotropic part of the brightness function during the tight coupling regime: 

• the radiation pressure given by the term /c^c^Gq which dominates on subhorizon 
scales and supports oscillations of the photon-baryon plasma; 

• the gravitational infall of the photon-baryon fluctuations in the potential well, A;^^'; 

• the ^> term that sources the oscillations at scales that enter in the horizon between 
the epochs of matter-radiation equality and decoupling. The $ term contribute as 
a friction term. 

Neglecting time variations in the gravitational potentials, Eq. 1)6. 22() describes a forced 
harmonic oscillator equation. On subhorizon scales the photon pressure balances the 
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gravitational collapse of the perturbations and causes the propagation of pressure waves in 
the photon-baryon plasma. The characteristic scale upon which these waves can propagate 
is fixed by the size of the sound horizon Vs = J Cs{rj)dri. When the photons last scatter, 
the imprint of these compression and rarefaction regions will appear as a series of peaks 
in the anisotropy power spectrum on the angular scales smaller than the angle subtended 
by the sound horizon at the decoupling. The odd and the even peaks correspond to 
photons coming from regions that are respectively in compression and expansion phases 
at last scattering. However the compression and rarefaction phases are not completely 
equivalent. In fact the baryons increase the effective mass of the barotropic fluid and 
consequently the gravitational infall generated by the density fluctuations leads to a deeper 
trough during the compression phase (baryon drag). The compression is enhanced over 
rarefaction resulting in a larger amplitude of the odd peaks relative to the even ones. 
Therefore the relative heights of the acoustic peaks is extremely sensitive to the baryon 
density Vtbh?. An additional source of anisotropics is due to the motion of the fluid at the 
last scattering surface due to the Doppler shift of the CMB photons. However due to the 
large value of R this effect is subdominant. 

Sachs- Wolfe effect 

In the large wavelength limit {krj << 1) the Boltzmann equation Eq. I|(i.l4|) reduces to 
the form Qq = — <I> ^I* (<I) = — ^f). The solution of this equation depends on the initial 
condition of the solution of the linear density perturbation equations. In fact Bq = 6r 
and for adiabatic initial conditions (^r(O) = — ^^'(0), reflecting the fact that the photons 
are overdensities inside the potential well. However due to the decay of the gravitational 
potential after matter-radiation equality the photon temperature at decoupling becomes 
0o(^*) = ('?*)) where ry* is the conformal time at decouplings. When the photons 
last scatter they climb out of the potential well ^ so that the effective superhorizon 
perturbations in the photon temperature is: 

[eo + *](r?.) = ^M'(7?.). (6.24) 

This is the so called Sachs- Wolfe effect [168]. It is the dominant source of anisotropies 
on large angular scales and is responsible for the plateau at the low multipoles in the 
CMB power spectrum. Since these scales are superhorizon at the time of decoupling 
the amplitude of this Sachs- Wolfe depends on the primordial spectrum of the density 
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fluctuations. 

The Integrated Sachs- Wolfe effect 

The time evohition of the gravitational potentials is a source of primary anisotropics, in 
fact as the photons free stream after decoupling the solution to the Boltzmann equation 
is given by [169]: 



where ji{kri) are the Bessel functions. The first term represents the contribution of the 
Sachs- Wolfe effect, while the second term is the Integrated Sachs- Wolfe effect (ISW) [170]. 
We can distinguish two contributions to the ISW: an early and a late ISW. The former 
is due to the decays of the gravitational potentials after horizon crossing at the end of 
the radiation dominated era. The latter occurs in dark energy dominated cosmologies, 
as the Universe starts accelerating the rapid expansion causes the decay of the density 
fiuctuations that drive the decay of the gravitational potentials. The combined effect of 
the SW and the two components of the ISW can be studied by analysing the the anisotropy 
power spectrum per logarithmic k and / interval defined by 



In figure IHTTl we show the density plots of Eq. (|6.26j) integrated up to z = 5 {left panel) 
and to z = {right panel) for a ACDM with = 0.7. The brighter regions correspond 
to values of the wavenumbers and multipoles where there is more power. Since for the 
model we have considered the acceleration starts at redshift z = 0.7, the contribution of 
the late ISW can be emphasized by comparing the distribution of power in the left panel 
with that in the right panel. As we expect, we note that due to the projection onto the 
last scattering surface the wavenumbers k and the multipoles I are strongly correlated. 
As a consequence of this the power at redshift 2: = 5 is shifted toward lower multipoles 
than at z = 0. This is simply because the last scattering surface is closer to the observer 
at redshift z = 5 and different scales are projected on larger angles. The filamentary and 
periodic structure of the density plots are due to the Bessel functions. We therefore expect 
a different structure in the case of an open or closed space-time. Note that the power is 
distributed along an upper and a lower ridge respectively. The latter is the contribution 
of the Sachs- Wolfe effect caused by modes that are superhorizon at decoupling {k ~ 10"'^) 




(6.25) 




(6.26) 
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Figure 6.1: Density plots of Ci{k) in tlie k — I plane for a ACDM model. The modes in the 
left panel are integrated up to redshift z = 5 and those in the right panel are integrated 
up to z = 0. 

while the upper ridge is produced by the early ISW. Since this latter effect arises closer 
to the present time, the corresponding anisotropics are projected on angular scales larger 
than the SW effect itself. The bright spot at the top of the panels correspond to the rise 
of the first Doppler peak, the other peaks do not appear in these plots since we limited 
our analysis to scales larger than the horizon at decoupling {k < 10~^). As we can see in 
the right panel, integrating the modes over the period of time when the Universe started 
to accelerate boosts the power at lower multipoles. In fact the signature of the late ISW 
can be noticed from the fact that the upper ridge becomes brighter than that in the right 
panel. Note that the two plots have slightly different brightness scales, this is because the 
late ISW overlaps with the SW leading to a different distribution of the power on the large 
angular scales. 

Damping mechanisms 

The overall amplitude of the CMB power spectrum is suppressed by the effects of diffusive 
mechanisms. The most important of these effects is due to the diffusion of the photons 
on scales smaller than their mean free path Ac*. In fact due to the Compton scattering 
the photons randomly walk through the baryons with a mean free path Ac ~ t^^- After 
N collisions the diffusive length is A/) ~ ^/NXc- Consequently on scales A < A/j the 
diffusion exponentially suppresses the amplitude of anisotropics. The baryons suffer a 
similar effect, in fact from Eq. (|6.12|) we can see that on scales k « f/R the Compton 
scattering can drag the baryons in and out of the potential wells leading to a destruction 
of the baryonic oscillations (Silk damping [171]). Another damping effect arises from the 
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fact that he last scattering process is not instantaneous. The photons we observe today 
in a given direction of the sky may have decoupled from different points along the line 
of sight. Therefore fluctuations on scales smaller then the thickness of the last scattering 
surface have destructively interfered causing a 'washing out' of the anisotropics on very 
small angular scales. The reionization history of the Universe also contributes to the over- 
all damping of the temperature anisotropy power spectrum. In fact when the first stars 
form they reionize the intergalactic medium, therefore a fraction of the CMB photons are 
scattered on scales smaller then size of the horizon at the epoch of reionization. This 
lead to a suppression of the power at large multipoles. In light of the recent WMAP 
measurements of the temperature-polarization cross power spectrum, it is worth men- 
tioning that the reionization sources the polarization of the anisotropies. Consequently 
a late reionization, Zion < 100, will produce significant power at low multipoles in the 
temperature-polarization spectrum. 



6.2 Dark energy and the Integrated Sachs- Wolfe effect 

We have previously discussed the importance of distinguishing between different dark 
energy models from the "concordance" ACDM model. The ISW is particularly sensitive 
to the late time evolution of the Universe. In fact during the matter dominated era the 
gravitational potential $ associated with the density perturbations is constant and there is 
no ISW effect. However, we have seen that in ACDM models $ starts decaying at redshifts 
when A starts to dominate, producing large angular scale anisotropies [172]. In dark energy 
scenarios the cosmic acceleration is not the only contribution to the decay of <I>: on large 
scales the clustering properties alter the growth rate of matter perturbations [76,173]. It 
is the signal of this clustering [174] that we are hunting in as model-independent a way 
as possible. We assume a flat spatial geometry and fix the value of the Hubble constant 
Ho = 70 Km.s~^ Mpc~^ , and the amount of matter (CDM) = 0.3. We can usefully 
distinguish two classes of models, 

• those with a slowly varying equation of state for which < aj? /A < 1 , as in the case 
of the inverse power law potential [79]; 

• a rapidly varying w{z), such as the 'Albrecht-Skordis' model [118] and the two expo- 
nential potential [100], with a™/ A > 1. This class also includes many interesting rad- 



6.2 Dark energy and the Integrated Sachs- Wolfe effect 



70 



0.2 




10° 10"' 10"^ 10"^ 

a 



Figure 6.2: Time evolution of the equation of state for two classes of models, with slow 
(red solid line) and rapid transition (blue dotted line). The dark energy parameters specify 
the features of wqi^a) 

ical models such as vacuum metamorphosis [175], late-time phase transitions [176], 
and backreaction- induced acceleration [177]. 

We show these two classes in figure ESI The red solid line corresponds to dark energy that 
tracks the dust during the matter era {wq = 0.0) and evolves slowly toward Wq = — 1, and 
the blue dotted line corresponds to a model with a rapid transition in its equation of state 
at a"^ = 0.1 (z = 9). Given current data it is worth studying the case with Wq = — 1, (since 
it is also the most difficult to distinguish from ACDM) whilst allowing the other parameters 
Wq and a™, to vary. We have modified the CMBFAST code [178] in order to include the 
effects of a dark energy fiuid and its perturbations by specifying the time evolution of 
the equation of state. Figure l6?3l shows the anisotropy power spectrum, C[^^ , produced 
through the integrated Sachs- Wolfe effect by a rapidly evolving {top panels) and a slowly 
evolving (bottom panels) equation of state; the red (solid) line corresponds to the ACDM 
model. As we can see in the top left panel (figure l6?3k ). varying a™ can produce a strong 
ISW. The effect is larger if the transition in the equation of state occurs at redshifts z < 3. 
On the other hand the Cf^^ is the same as in the cosmological constant regime if a"^ < 0.2 
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Figure 6.3: Power spectrum of the ISW for rapidly varying models (top panels) and 
slowly varying ones (bottom panels). The solid red line is the ISW effects produced in the 
cosmological constant case. Detailed explanation in the text. 

(z > 4). In the top right panel ffigure l^^^ ) we plot the ISW for two different values of Wq, 
corresponding to Wq = 0.0 (dashed line) and Wq = —0.1 (dot-dash line). We note that the 
signal is larger if the quintessence field is perfectly tracking the background component. 
But as Wq diverges from the dust value the ISW effect approaches that of ACDM. This 
means that even for rapidly varying w{z) (small A), the ISW is distinguishable from that 
in the ACDM scenario only if w{z) during matter domination closely mimics the dust 
value and the transition occurs at low redshifts, z < 3. We can see that the amplitude of 
the integrated Sachs- Wolfe effect is smaller in slowly varying models (bottom panels). As 
we expect the Cj^^ is independent of a™ (figure EIHt), since for these models a different 
value in the transition redshift does not produce a large effect on the evolution of the 
dark energy density. In figure IH3b . the ISW power spectrum is large for Wq = 0.0 (dash 
line) and becomes smaller than the cosmological constant on horizons scales as Wq has 
negative values (dot-dashed line), and increases toward A for Wq approaching —1. This 
class of models is then more difficult to distinguish from the ACDM if the equation of 
state today is close to = — 1- This can be qualitatively explained noting that perfect 
tracking between dark energy and CDM causes a delay in the time when the gravitational 
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potential starts to decay, compared to the case of ACDM. This effect is stronger for models 
with a rapidly varying equation of state since the rapid change in wq produces a stronger 
variation in the gravitational potential. 



6.3 Differentiating dark energy models with CMB measure- 
ments 

The information in the CMB power spectrum can be encoded in the position of the first 
three Doppler peaks, li, I2 and ^3 and in three parameters. Hi, H2 and that account 
for the amplitude [179]. These are the height of the first peak relative to the power at 
/ = 10, 

the height of the second peak relative to the first, 
and the height of the third peak relative to the first, 

In principle the position of the second and third peaks is not necessary to characterize the 
CMB power spectrum, since their value is set by the position of the first peak through the 
harmonic relation Eq. H4.6() . However, as we have previously discussed in Chapter |1] this 
relation is affected by pre-recombination effects and therefore the value of I2 and ^3 carry 
information about a combination of the cosmological parameters. Since Hi, H2 and 
quantify the amplitude of the power spectrum on different multipoles we expect they are 
sensitive in different ways to the cosmological parameters. For instance since Hi depends 
on the ratio between C/^ and C/^^ it mainly depends on the scalar spectral index n, the 
physical baryon density O^/i^ due to the baryon drag that changes Cj^, Qa through the 
ISW and the reionization since increasing r will lower the first peak. On the other hand H2 
is sensitive only to Q^h'^ and n. Since the baryons affect in the same way the height of the 
third and first peak, ^^3 depends only on the matter density and n. This implies that 
if the amplitude of the first three peaks is accurately determined the degeneracy between 
^l^h"^ and n can be broken and their value can be measured with high precision. We use 
Hi, H2 and H3 to determine the overall effect of the dark energy in the CMB spectrum. 
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Figure 6.4: Relative difference of H\ (blue solid line), (green dash line) and ifa (red 
dash-dot line) to the ACDM model, for rapidly varying models (top panels), and with 
slow transition (bottom panels). For these models the present value of the equation of 
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In particular we expect Hi to be extremely sensitive to a quintessential effect through the 
ISW, while H2 and H3 are weak indicators of such a signal. We have computed the CMB 
spectra for the class of models with Wq = — 1 described in the previous section and we 
have inferred the values of Hi and compared them with those of the ACDM model for the 
same values of the cosmological parameters. In figure l6^ we plot the absolute value of the 
relative difference of Hi, H2 and H3 to the ACDM model. The rapidly varying models are 
shown in the top panels. We can see the strong ISW effects produced by changing a"^ are 
now evident in the large discrepancy between Hi and H^ (blue line) (fierure lH^lK ): it can be 
larger then 20 per cent for a™ > 0.6. The effect on H2 and H-^ is smaller. However varying 
Wq (figure 16.4b ') produces a discrepancy of only order 4 per cent on Hi , while H2 and 
H^ remain the same as in ACDM. For a slowly varying equation of state. Hi, H2 and H3 
are independent of a™ (figure EUt) . The dark energy imprint is only on Hi for which the 
discrepancy to the A case is about 10 per cent. This discrepancy decreases when changing 
the value of the quintessence equation of state during matter from Wq = 0.0 to Wq = — 1 
f figure l6^1h ) . Values of the equation of state today Wq > — 1 imply a stronger ISW effect. 
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Figure 6.5: As in figure HTH for Wq = —0.88. 

Consequently tlie curves of fig. 3 are shifted upwards. For instance in figure 16.51 we plot 
the class of models previously analysed, with Wq = —0.88. We note the same behavior 
as we vary the dark energy parameters, but the discrepancy with the ACDM model is 
now larger. In figure IH3b it is worth noticing the case = — 1, that corresponds to a 
model very similar to a 'k-essence' model [157]. We can see that the relative difference 
with the ACDM case is of the order of a few percent, in agreement with [180] for the same 
value of Wq = —0.88. At this point we ask the key question whether such differences are 
observable. We have shown that Hi is a good estimator of the ISW effect, and that it is a 
tracer of the dark energy imprint on the CMB. However its estimation from the data will 
be affected by cosmic variance at / = 10. Hence with even perfect measurements of the 
first acoustic peak the uncertainty on Hi will be dominated by the 30 per cent uncertainty 
due to cosmic variance. With the plots of figures I?mi6.5l in mind, this means that if the 
present value of the equation of state is close to —1, slowly varying dark energy models are 
hardly distinguishable from ACDM, while rapidly varying ones can produce a detectable 
signature only if the transition in the equation occurred at a™ > 0.7, but in any case it 
will be difficult to constrain Wq. Moreover it should be taken into account that Hi is 
degenerate with other cosmological parameters, such as the scalar spectral index n, the 
optical depth r and the scalar to tensor ratio r. Hence only an accurate determination 
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Figure 6.6: Relative difference of l\ (blue solid line), I2 (green dashed line) and Z3 (red 
dot-dashed line) to the ACDM model, for rapidly varying models (top panels), and with 
slow transition (bottom panels). For these models the present value of the equation of 
state is Wq = — 1. 

of the angular diameter distance, inferred from the location of the acoustic peaks, would 
allow detection of such deviations from the cosmological constant model. The shift of the 
multipole positions (ii) of the acoustic peaks caused by the evolution of the dark energy 
in the class of models analysed in figure 16.31 can be seen in figure 16.61 where we plot 
the relative difference of li,l2 and Z3 to the A case. We note that due to the additional 
shift induced on the first acoustic peak by the ISW effect the difference with the ACDM 
model for the first peak is generally larger than for the second and third peaks. As with 
the comparison of the amplitude of the CMB spectrum, the largest effect is produced by 
models with a rapid transition occurring at small redshifts. However the degeneracy of 
the angular diameter distance, in particular with the value of Hubble constant and the 
amount of dark energy density, will limit our ability to put tight constraints on the dark 
energy parameters. There are alternative ways in which these problems can be alleviated, 
for instance cross-correlating the ISW effect with the large scale structure of the local 
universe [181-183]. An efficient approach would be to combine different observations in 
order to break the degeneracies with the cosmological parameters [184,185]. 
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Figure 6.7: CMB power spectrum for a fiducial dark energy models, the errorbars are the 
cosmic variance errors. 

6.4 Testing dark energy with ideal CMB experiments 

The new generation of CMB satellite experiments is going to provide an estimation of the 
anisotropy power spectrum close to the ideal case, where the only source of indetermination 
is due to the cosmic variance. It is therefore interesting to test the sensitivity of such CMB 
measurements to the dark energy effects. As we have described in the previous section, the 
class of models most difficult to distinguish from the ACDM case corresponds to a dark 
energy fluid with a slowly varying equation of state characterized by Wq = — 1. This can 
be considered as the most pessimistic situation since it can prevent us from understanding 
the nature of the dark energy. We have studied the information on this class of models that 
can be inferred from cosmic variance limited measurements [162] and generated a sample 
of ideal CMB power spectrum data assuming a fiducial model specified by the following 
values of the cosmological parameters: Q,q = 0.68, Hq = 70 Km Mpc~^, Wq = — 1, 
Wq = —0.4, = 0.2, A = 0.3. Assuming a fiat geometry, no tensor contribution, 
the scalar spectral index n = 1 and the baryon density O^/i^ = 0.021, we have binned the 
input power spectrum in 62 data points plotted in figure IHTl A library of CMB spectra has 
then been generated using a modified version of CMBFAST [178] assuming the following 
uniform priors: G (0.5,0.8), Hq G (0.5,0.8), w% G (-0.1,-0.4), u;™ G (-1.0,0.0) and 
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Figure 6.8: Simulated 2-D likelihood contour plots for an ideal experiment, with no priors 
{top panels) and Hq prior (bottom panels). The yellow and red contours correspond to 1 
and 2(7 respectively. 

a™ G (0.01, 1.0). As extra priors we have assumed A = 0.5 (slowly varying equation of 
state) and the remaining parameters have been set as follows: n = 1, = 0.05 and r = 0. 
Since we are not considering the effect of systematics, but only cosmic variance errors we 
evaluate a simple likelihood defined as 

Cf{k;aj)-Cf{k 



L{aj) = jV exp 



(6.30) 



where aj are the likelihood parameters, jV is a normalization constant and o", is the 
cosmic variance at I = li. The superscripts t and d refer to the theoretical quantity and 
to the real data respectively. The results are shown in figure ESI where we plot the two 
dimensional likelihood contours in the Wq — CIq (left panels) and Wq — a™ (right panels) 
planes respectively. We find Qq = 0.68ibQQg, Wq < —0.78, whereas assuming a prior on 
Hq improves the constraints on Wq to Wq < —0.85 at la. The likelihood plot in the 
Wq — a™ plane shows that is undetermined. This result is expected since for this class 
of models the value of a™ does not affect the evolution of the dark energy. On the contrary 
we find w"^ = -0.79±[|:f, but the A case cannot be excluded at 2a. The constraint does 
not improve assuming the Hq prior. This is because Wq and Wq are degenerate, hence 
marginalizing the likelihood over Wq shifts the best fit value of Wq towards more negative 
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values. However it is remarkable that there is still some sensitivity to the value of Wq, 
such that the best fit is not for Wq = — 1. Therefore we can conclude that if Wq ~ — 1, 
a large class of models will not be distinguished from a ACDM scenario even with ideal 
CMB measurements. It is possible that by combining different cosmological data, as Sn 
la, large scale structure and quasar clustering the degeneracy between Wq and Wq can be 
broken and more information on this class of dark energy models can be inferred [185]. 



Chapter 7 

Alternative cosmological test with 

higher order statistics 



In the recent past the number of papers devoted to the analysis of high order statistics of 
the CMB anisotropy has dramaticahy increased. In fact the simplest inflationary models 
predict to first order a Gaussian distribution of temperature fluctuations, and deviations 
from gaussianity could be the signature of other phenomena occurring in the anisotropy 
formation process. On top of that different mechanisms can also be a source of non- 
Gaussian signals at different angular scales. This is the case of non-linear effects during the 
inflationary epoch or the presence of topological defects. The present CMB data strongly 
constrains the level of non-gaussianity on the scales so far probed by the experiments. In 
this Chapter we will briefly review the higher order statistics of the CMB anisotropics. 
We will focus on the use of the angular bispectrum as an estimator of non-gaussianity, and 
introduce a formalism that allows us to analytically calculate the spectrum and bispectrum 
in the case of a random distribution of localized anisotropics in the CMB sky. We will argue 
that applying this analytical approach to the analysis of localized anisotropics such as the 
Sunyaev-Zeldovich effect or radio point sources, it is possible to constrain the clustering 
properties of these objects and determine the cosmological parameters. 
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7.1 Higher order statistics 

Let us expand the temperature fluctuation field in the direction 7 of the sky into spherical 
harmonics: 

rj. =2^ aimYi (7), (7.1) 

l,m 

the aim are the multipole coefficients that contain all the statistical information of the 
anisotropy field. For m ^ these are complex numbers satisfying the condition a^-m = 
(— l™')a^^. The statistics of the CMB anisotropies depends on the physical process re- 
sponsible for the generation of the initial density perturbations. Simple infiationary mod- 
els predicts a Gaussian spectrum of fiuctuations [186, 187], consequently we expect the 
anisotropy field to be Gaussian. In such a case the aim are random Gaussian variables 
with Gaussian distributed amplitudes and with uniformly distributed phases. As conse- 
quence of this the statistical distribution of the CMB anisotropies is entirely specified by 
its second order moment (i.e. the power spectrum), 
I 

m=—l 

In this Gaussian case the higher odd moments of the distribution vanish, while the even 
moments can be expressed in term of the variance Ci. Therefore any deviation from 
gaussianity will inevitably produce non vanishing high order statistics. Since there is a 
potentially infinite number of higher moments, the non-Gaussian hypothesis cannot be 
disproved. Several methods have been proposed in the literature to measure statistical 
estimators of the skewness (third moment) and (kurtosis) from the analysis of CMB maps 
(for a list of these methods we refer to [188]). A simple method is to measure the correlation 
between the temperature fiuctuations in different directions of the sky. In this case the 
higher order moments are estimated by the angular correlation functions. For instance 
the m — point angular correlation function is defined as: 

/AT AT AT \ 

Cm(7l,72,-",7m) = ( — (71) — (72)-"^(7,n) ) , (7.3) 

where 71, .., 7m are unitary vectors pointing at 1, .., m directions of the sky and the average 
(..) is taken over the whole sky. In particular an estimate of the skewness is provided by 
the collapsed three-point correlation functions C^{a), that is a specific configuration of 
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the three-point correlation function between two points of the sky separated by an angle 

a and defined as: 

/AT AT AT \ 
C3(a) ^ ^ — (^0 — (71)^(72)^ , (7.4) 

where cos(q;) = 71 • 72- In terms of the multipoles it reads as: 

^3(a)= E ^/l(cosa)a^,^,a^,„,a;3^3#^l#^2#^3J^™"^^ (7.5) 

where P;j(cosq;) is the Legendre polynomial of degree li, Wi is the experimental window 
function in the multipole space and 

^m,m,m, _ ^ / (2£l + l)(2£2 + l)(24 + l) ,3 X 

are a combination of the Wigner-3J symbols. The C^{a) measured from the COBE- 
DMR sky maps has been found to be consistent with the Gaussian hypothesis within 
the cosmic variance errors [188]. These results strongly limit the allowed amount of non- 
gaussianity in the CMB anisotropies at large angular scales, in particular they constrain the 
parameter space of a class of non-Gaussian models [189-192]. Similarly the collapsed three- 
point correlation function inferred from the WMAP data is also consistent with Gaussian 
expectations [193]. Any non-gaussian analysis can be carried out in the multipole space 
as well, in this case the equivalent of the three-point correlation function is the angular 
bispectrum defined by: 

^hhh — 2-^ \m.x m.2 m-i J "'l^ % "«3 • \'-' ) 

mi ,m2,m3 

This estimator is rotationally invariant [194, 195] and satisfies geometrical conditions such 
that the only non vanishing Bi^i^i^ are those with |/j — lj\ < Ik < h + Ij (i 7^ i 7^ k) for 
all permutations of indices and h + h + ^3 = even. A normalized version of the angular 
bispectrum has been applied to the analysis of COBE-DMR maps for the configuration 
h = h = ^3 [196-199]. More recently the bispectrum has been estimated from the WMAP 
map and has been found to be consistent with the Gaussian hypothesis [200]. 

7.2 Frequentist approach and estimation of higher moments 

The frequentist approach is the usual procedure adopted to test gaussianity. This im- 
plies the measured value of a higher order statistical estimator is compared against the 
probability distribution function obtained from random Gaussian simulations of the data 
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sample. If there is a low probability that the measured value is consistent with the Gaus- 
sian simulation, then the Gaussian hypothesis is ruled out. For instance this approach 
allowed the authors of [196] to rule out the gaussianity of the anisotropies at multipole 
Z = 16 in the COBE-DMR data. However it is worth mentioning that so far only the 
diagonal component of the bispectrum has been measured and in order to be statistically 
significant the data analysis should be extended to the estimation of the non-diagonal term 
of the bispectrum. If this is not a problem at low multipoles it could be a computationaly 
challenge at higher orders. On the other hand it has been pointed out in a number of 
papers [199, 201] that within the frequentist approach some non-Gaussian theories will be 
indistinguishable from the Gaussian one. We will try to make this point more clear with 
a specific example and we refer to the fundamental statistics textbook [202] for a more 
detailed derivation of the formula used in what follows. 

Let us consider a sample of data {xj} (i = 1, .., N) generated from a random Gaussian 
process f{x) such that each Xi is an independent Gaussian random variable with zero 
mean /ii = and the second order moment fi2 = 2, 

/(x) = -^e-^. (7.8) 



In order to avoid confusion we will use Greek letters for the moments of the generating 
function /(x), fir its r-th moment and Roman letters for the moments obtained from the 
statistics of the data sample. In particular the r-th moment rrir is given by: 
1 " 

mr = j^J2('^i-{x)Y, (7.9) 

i=l 

where (x) = mi is the mean value of the sample. Making K realizations of this data 
sample, {xi}j (j = 1, .., K), we can infer the distribution function Pgauss{fnr) of the r — th 
moment statistic m^, which is the frequency (number of times to the total number of 
realizations) of a certain value appearing in the sample Xij. The expectation value of 
this distribution E{mr) (the mean value of estimated from the K realizations) will be 
indicative of the r-th moment fi^ of the generating Gaussian distribution function. The 
expectation value E[mj.) and the variance var{mr) can be related to the moments /x^ of 
the generating function through approximate relations valid up to (K ■ N)^^/"^ [202]. It 
can be shown that: 

E{mr) = fir- (7.10) 

var{mr) = -r^—rrXl^2r - f^l + r'^fi2f4-i - 2r^r_i/ir+i). (7.11) 
iV • K 
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For practical purposes we prefer to work in terms of cumulants kj., that are an equiva- 
lent set of numbers characterizing the generating function f{x). The cumulants kr and 
the moments fir ai'e related through Eq. (3.33) in [202]. Let us consider the third or- 
der moment-statistic m^, from Eq. 1)7. 10(1 and Eq. H7.11() and using the relation between 
moments and cumulants we have: 

E{ms) = ks, (7.12) 
varims) = J^ih + Qk2k^ + Gkl + 9A:|). (7.13) 

For the Gaussian generating function we have considered A;2 = /-f2 = k^ = 0, k^ = 

and kQ = 45(7^, consequently we obtain: 

E{m3) = 0, (7.14) 
183 

var{m3) = 4^.^ ^^- (7-15) 

On the other hand let us consider a sample of data generated from a non-Gaussian random 
process. A general non-Gaussian random generating function can be constructed by using 
an Edgeworth expansion around the Gaussian distribution defined by 

g{x) = f{x){l + hsHsix) + ^k^H^ix) + ...), (7.16) 
D 24 

where f{x) is given by Eq. 1)7. 8|) . the Hi are the Hermite polynomials and the cumulants 
kr are free parameters. It can be shown that Eq. ()7.16)) gives a good approximation to 
any distribution function provided all moments are defined and the higher order terms do 
not dominate over the Gaussian one. We can limit the non-gaussianity to first order by 
imposing only the presence of a non vanishing skewness k^ and kurtosis /C4. In this case 
Eq. (TrTUI) and Eq. ^TT^i becomes: 

E{ms) = ks, (7.17) 

var{ms) = ^(^^4^' + + 9^1)- (7-18) 

Note from Eq. 1)7.18(1 the kurtosis and the skewness can sum up in such a way that 
they cancel each other. In such a case the frequency distributions of the third order 
moment-statistic m-s inferred from the K realizations of the Gaussian and non-Gaussian 
process, Pgaussi'^^'i) a-iid Pn-gauss{''^z) respectively, will have the same variance. This 
makes difficult to establish whether the population of the data had been generated from a 
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Figure 7.1: Probability distribution functions of tlie third order moment estimated 
from Monte-Carlo simulation of randomly generated data for a Gaussian (red bars) and a 
non-Gaussian process (blue bars). 

Gaussian or a non-Gaussian process by estimating ms. This analysis is confirmed by the 
numerical simulations in which we simulated two populations of data from two different 
random processes. The Gaussian sample has been generated assuming the generating 
function Eq. (|7.8|) with ki = 0, /c2 = 1, while the non-Gaussian population has been 
obtained assuming Eq. (frT6|) with ki = 0, k2 = I, h = 3 and k^ = -153/18. From these 
two populations we inferred the distribution functions of the third moment-statistic ms. 
The results are shown in figure ITTI where we plot Pgaussi^s) (fed bars) and Pn-gaussi^s) 
{blue bars). Note that Pn- gauss {ms) is peaked at ma ^ 0, while Pgaussim-s) has a maximum 
at ms = 0. Even so the variances of the two distributions are the same. This naive 
argument shows that for some non-Gaussian models little information can be obtained 
using a frequentist approach when analysing the third order statistics unless higher order 
moments of the CMB anisotropies fields are measured as well. Such simulations are 
necessary to take into account the systematic experimental sources of non-gaussianity. 
However they make these tests computationally expensive. Therefore it is important to 
have a theoretical prediction of the amplitude of non-Gaussianity at least for known non- 
Gaussian anisotropies such as the secondary ones. 
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7.3 Modelling localized non-Gaussian anisotropies 

Localized patterns of anisotropies produce a non-Gaussian signature in the CMB statistics. 
This is the case of secondary sources of anisotropies such as radio point sources or the 
Sunyaev-Zeldovich (SZ) effect caused by the interaction of the CMB photons with the hot 
gas associated with clusters of galaxies [203]. Active sources of non-Gaussianity include 
topological defects such as cosmic strings, domain walls (see for a general review [204]) or 
primordial bubble relics of a first order phase transition during the inflationary era [205, 
206] . For these classes of models analytical formula of the higher order correlation functions 
have been calculated in a number of papers [191, 192,208]. As far as the non-Gaussianity 
arising from secondary source of anisotropies is concerned, an analytical estimate of the 
three point correlation function has been obtained in the case of the Rees-Sciama effect 
[207], while the value of the bispectrum produced by the S-Z effect and extragalactic 
radio sources has been obtained in [209], whilst the case of the Vishniac effect has been 
considered in [210]. In what follows we determine the spectrum and the bispectrum for 
the case of a distribution of localized anisotropies in the CMB sky, using a formalism 
introduced in [208]. 

Power Spectrum 

Let us consider a nearly circular spot on the sky and perform the spherical harmonic 
decomposition in the frame where the z-axis coincides with the center of the spot. The 
brightness of the temperature fluctuation of the spot is b, and f{9) is its angular profile. 
In this frame, we can write the temperature fluctuation as As{0) = b ■ f{9), expanding it 
in spherical harmonics we find [208]: 



where Pi{cos9) is the Legendre polynomial and eim is a perturbation induced by the irregu- 
larity of the spot. As a first approximation it can be neglected, as the CMB measurements 
do not have the enough resolution to resolve its (/>-angular structure. For a distribution 
of N spots, performing a rotation to a general frame where the n-th spot is at the angle 




(7.19) 



with 




(7.20) 
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In = {dn-, 4>n) and summing over all the spots we obtain: 

aim = Y.^n^nyr^in), (7.21) 

n 

where 6„ is the brightness of the n-th spot and 

yi = 2TT J d{cOSe)fn{e)Pl{cOSe) (7.22) 

is the Legendre transform of its temperature profile. We can now calculate the power 
spectrum, 
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n n' 



21 + 1 



= E E ^nV^n^n'C^nn' " (7-23) 
n n' 

where we have used the completeness relation of the spherical harmonics and Cnn' = 
Snn' + Wnn', with Wnn' taking into account the possibility that the spots are uncorrelated 
{6nn'), as for a Poisson distribution, or that are correlated two by two {wnn')- Therefore 
by expanding the sums we finally obtain: 



^bliJ^nf + ^ E bnbn'Wnn'^ij^i'Pliinfn'] 



2 



(7.24) 



We assume now that the spots have about the same brightness and the same temperature 
profile: b = (bn) and / = (/„) where the average is taken over the ensemble of spots. 
Because the number of spots in a circular ring centered on a single spot is proportional 
to the angular extension of the ring we can substitute the sums over the number of spots 
with an integral over the whole sky. Hence we can substitute the discrete correlation 
coefficient Wnn' with the two-point angular correlation of the spots, w{a). If the spots 
correspond to radio sources or clusters of galaxies w{a) would be the angular correlation 
function estimated from large scale structure observations. The mean value over the spot 
distribution of Eq. 1)7. 24() becomes: 

{Ci) = Nb^ + Ji ■ ^i] , (7.25) 

where is a normalization constant and 

'^i = J d{cosa)J^i{a)wia)Pi{cosa), (7.26) 
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with 

^i{a)= j d{cose)f{9 + a) Pi {cos 9). (7.27) 

From Eq. (|7.25|) we may note that the contribution to the anisotropy power spectrum 
from a localized distribution of spots increases with the brightness b. If the spots are 
uncorrelated the second term J^i ■ drops from the equation and the only contribution 
is due to the projection of the spot's signal in the multipole space J^i. In particular the 
largest contribution will occur in the range of multipoles that correspond to the effective 
angular size of the spots. If the spots are correlated, the overall contribution to the power 
spectrum will depend on the sign of J'l-'^i. 

Bispectrum 

Following the same procedure as just described we calculate the angular bispectrum. Sub- 
stituting Eq. (frnj) in Eq. ifTT)) we have: 



-^hhh — 2-^ \mx m2 ma j "ii h "«3 

mi ,m2,m3 

!,™3 

E ft/4'4)>^r*(7nj^r*(7n2)iir(^"3), (7.28) 



ni,n2,n3 
X 



mi,m2,m3 

where Cni^n^no, takes into account all the possible correlations between the spots up to 
third order, 

Cn-yn^n-i — ^n\n2^n-2,n-i ~l~ '^nxn^^n^n^ ~l~ ^nin3'^n2n2 ~l~ '^n^n-i^nxrix ~l~ '^nxn^n-i,^ (7.29) 

with Wnin2n3 the correlation between three different spots. Using the definition of the 
Wigner 3J symbols in terms of an angular integral of three spherical harmonics and the 
completeness relation of the spherical harmonics Eq. ()7.28|) becomes: 



o 

Eh3 ah ah ah . a2>hhh j u ah ah ah^,, 
"ni=^ni^ni=^ni ^^7^^ ^ 2 ^ "ni ""2 =^ni =^ni =^n2 "^"li^a 

"■1 ni^n2 



^h ah ah.,, . aphhh 



(7.30) 



where 



V(2Zi + l)(2/2 + l)(2/3 + l) ^° 00/' ^ ^ 



7.4 Discussion 



88 



csihhh 



dn^Pl^ (7 • 7ni )^«2 (7 • lni)Pl3 (7 -7^1), 



(7.32) 

(7.33) 
(7.34) 



Due to the background isotropy of the space, we are free to choose 7 = . As consequence 
the integrals Eq. (|7.32I7.34|) become: 



'^7ni7n2 

7^]^ 7^2 7^3 



47r, 



'^5^7ni^«3(7 •7n2), 

d^ln^ P12 (7 • 7n2 )^i3 (7 • 7n3 ) • 



(7.35) 
(7.36) 

(7.37) 



As for the power spectrum, we can average over the ensemble of all the spots and substitute 
the sums over the spots with integrals over the whole sky. In this case the discrete correla- 
tion coefficient if^n-nj and WmujUk are replaced with the corresponding angular correlation 
function w^a) and w{a,f3). After tedious calculations we obtain: 



^(0) + 3 (1) 1 (3) 

2 3 



where 



(7.38) 

(7.39) 
(7.40) 

d(cos /3) J^Z2 (")=^/3 iP)w{a, (3)Pi^ (cos a)Pi^ (cos /3) , (7.41) 

with cJ^i and ^^^3 defined by Eq. (|7.22j) and Eq. (|7.26j) respectively and =/;(«) defined by 
Eq. ()7.27p . As we can see from Eq. ()7.38p for a Poisson distribution ^^^^ and ^^^^ vanish, 
but the bispectrum remains non-vanishing due to the term I^^^^ that account for the 
localized structure of the anisotropics. 



^(2) = j^i^ f d{cos a) 



7.4 Discussion 

In the previous section we presented a general formalism to calculate the contribution to 
the power spectrum and the bispectrum of a distribution of spots in the CMB sky. The 
formulae Eq. (|7.25j) and Eq. (|7.38|) have to be considered as a starting point for further 
investigation. They take into account several effects, for instance the brightness and the 
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angular size of the signals contribute at zeroth order to the spectrum and the bispectrum. 
On the other hand the presence of internal correlations in the spot distribution, which are 
described by the angular correlation functions t(;(a) and w{a,l3), contribute as first and 
second order effects respectively. As a specific application, Eq. ()7.25|) and Eq. ()7.38() can 
be computed in the case of a distribution of spots caused by the SZ effect of a cluster of 
galaxies. Approximating the shape of this signal with a Gaussian profile characterized by a 
given width, it will be possible to numerically compute the integrals Eq. (|7.22() . Moreover 
without loss of generality the angular correlation function w{a) can be assumed to be a 
power law. In such a case the integrals Eq. (|7.26j) can also be numerically computed. The 
resulting power spectrum can be compared with the prediction of numerical simulations. 
This will allow us to test how crucial is the assumption that b = and / = (/«). 
The next generation of CMB measurements will measure the bispectrum at very high 
multipoles, therefore using Eq. 1)7. 38(1 in a specific case such as the SZ effect offers an 
alternative way of inferring cosmological information. 



Conclusion and prospects 



In this thesis we have discussed various aspects of dark energy dominated cosmologies. In 
Chapter 1 we have reviewed the observational evidence of the dark energy. We have seen 
that different cosmological measurements are consistent only if the dark energy accounts 
for most of the matter content of the Universe. In Chapter 2 we have discussed some of 
the proposed dark energy candidates and we have focused on the quintessence scenario. 
In spite of the theoretical difficulties of a viable quintessence model building, this scenario 
has a number of interesting features that can be tested with cosmological measurements. 
In particular in Chapter 3 we reviewed the dynamics of scalar field perturbations for two 
different class of minimally coupled quintessence models. We have learnt that quintessence 
perturbations have no active role during the structure formation. However their presence 
can lead to time integrated effects in the evolution of the gravitational potential. Therefore 
they can leave a characteristic imprint in the Cosmic Microwave Background anisotropy 
power spectrum through the Integrated Sachs- Wolfe effect. In Chapter 4 we have con- 
strained with Sn la data and the position of the CMB peaks a parameterized quintessence 
potential that accounts for a large class of quintessence models. Using the properties of 
the tracker regime we have been able to put upper limits on the present value of the 
quintessence equation of state. We have found that by the present time the scalar field 
is evolving in flat region or close to a minimum of its potential. However the results of 
this analysis clearly indicate that the possibility the quintessence equation of state was 
largely different from its present value cannot be excluded. In this direction a lot of effort 
has been made to detect time variation of the quintessence with cosmological distance 
measurements. In Chapter 5 we reviewed some of the proposed methods and we pointed 
out a number of potential problems. In particular the use of a constant equation of state 
to parameterize the dark energy leads to misleading conclusions. Using a very general 
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argument we showed that if the dark energy is time varying the constraints on a constant 
equation of state will be pushed towards large negative values. Therefore all the results 
obtained using this approach have to be carefully interpreted. In particular the fact that 
several data analysis found the equation of state best fit value to be ii; < — 1 can be just a 
bias effect. On the contrary we have proposed a time parameterization of the dark energy 
equation of state in terms of physical parameters. This accounts for most of the proposed 
dark energy models and moreover is valid at all the redshifts. Hence this approach allows 
us to take into account in a model independent way not only the effects dark energy has 
had on the expansion rate of the Universe but also on the structure formation. In Chap- 
ter 5 we have applied this parameterization to study the dark energy effects in the CMB 
power spectrum. As conjectured in Chapter 3 we have found that the dark energy leaves 
a characteristic imprint through the ISW effect. The amplitude of this signature selects 
the class of models which are distinguishable from the cosmological constant scenario. In 
particular we have shown that using ideal CMB data, only models characterized by a rapid 
transition of the equation of state can be distinguished from the A case. In Chapter 7 we 
have introduced an alternative cosmological test using higher order statistics of the CMB 
anisotropies. We have computed the spectrum and the bispectrum for a distribution of 
localized non-Gaussian anisotropies. These can be applied to specific cases such as the 
SZ effect to constrain cosmological parameters through the non-gaussianity produced by 
the imprint of cluster of galaxies. We can find a number of directions where the work so 
far reviewed can be further extended. Under some general assumptions it would be inter- 
esting to test the formulae developed in Chapter 7 with the predictions of the SZ effect 
from numerical simulations of cluster galaxies [211]. A complete likelihood analysis of the 
parameterized dark energy equation of state is currently in progress. We make use of the 
full cosmological data so far available. We intend to extend this analysis to the quasar 
clustering, that is a good candidate for testing the dark energy [212]. An interesting issue 
arises from the recent WMAP data. In fact it has been found that the quadrupole and 
octupole are suppressed in contrast with the prediction of ACDM cosmologies [32]. The 
possibility that this is caused by a cancellation mechanism between the SW and ISW 
effects due to clustering properties of the dark energy needs further investigation [35]. 
At moment there are no final conclusions about the nature of the dark energy, luckily 
the upcoming and future cosmological data will provide a new insight of the dark energy 
phenomenology. It has to be hoped that such measurements will help us to formulate 
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the new paradigm of Cosmology that wih ahows us to correctly address the dark energy 
problem. The history of science shows that no scientific activity is possible in subjects 
where no paradigms have been found. Therefore we should ask ourself what direction will 
the cosmological investigation take if the dark energy problem remains unsolved. 
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